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The kinetic Glauber-Ising spin chain is one of the very few exactly solvable models of non- 
equilibrium statistical mechanics. Nevertheless, existing solutions do not yield tractable expressions 
for two-time correlation and response functions of observables involving products of more than one 
or two spins. We use a new approach to solve explicitly the full hierarchy of differential equations 
for the correlation and response functions. From this general solution follow closed expressions for 
arbitrary multispin two-time correlation and response functions, for the case where the system is 
quenched from equilibrium at TS > to some arbitrary T > 0. By way of application, we give 
the results for two and four-spin two-time correlation and response functions. From the standard 
mapping, these also imply new exact results for two-time particle correlation and response functions 
in one-dimensional diffusion limited annihilation. 
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INTRODUCTION 



Glauber did some of his pioneering work in the field of nonequilibrium statistical mechanics when 
he introduced the kinetic version of the one-dimensional Ising model [jj] about 40 years ago. He 
made a very careful choice for the model dynamics which, while resembling characteristic features of 
ferromagnets, also admits an analytic approach. Amongst other things Glauber actually derived the 
general solutions for one and two-spin correlation functions in his original paper. 

Since then the Glauber-Ising model has continuously appeared in the literature in various contexts. 
For obvious reasons it has been studied in much detail as a toy model for ferromagnets ElBL H| . 
A link to one-dimensional diffusion-reaction processes has extended its relevance furtherp, |7], |Sj , and 
the model has also been employed as an abstract metaphor for tapped granular media The most 
recent occurrences in the literature relate nonequilibrium fluctuation-dissipation relations and related 
questions in higher-dimensional Ising models at criticality [ToL [Til IT^ | . 

A large fraction of the publications that deal with the Glauber-Ising chain either make direct use 
of the results derived by Glauber [2| , or use facts which can be obtained directly using his approach. 
This is rather remarkable, seeing as (at least) three different versions of general solutions for arbitrary 
order one-time spin correlations have been derived in the mean time. The authors of showed that a 
particular combination of correlation functions, denoted the C ra - functions, satisfy simpler dynamical 
equations than the correlations themselves. From the general solution given for the C„-functions, 
arbitrary (one-time) correlations could be obtained. As pointed out by the authors, however, their 
general expressions are rather complicated already for the case of three-spin correlations. In |l3| , the 
rigorous general solution of the master equation was obtained via a Fermion mapping. There, arbitrary 
correlation functions can in principle be found by expanding the associated operator according to a 
given scheme. However, as pointed out in the article, for correlations of more than two spins the 
amount of algebra involved increases very rapidly, and the formalism is not suitable for extracting two- 
time quantities. Finally, a rigorous and exact solution for the generating function of correlators has 
been found more recently using Grassmann variables. The formalism is very compact and yields an 
explicit expression for the generating function. Naturally, correlation functions are encoded implicitly 
in this result and follow by taking appropriate derivatives, but multispin or two-time quantities 
are again very difficult to obtain. The domain size distribution after quenching the system from a 
completely uncorrelated initial state to zero temperature is also briefly discussed in |l4j . The latter 
quantity was first obtained in |l5j using yet another method, namely the analogy between random 
walk problems and zero temperature dynamics of one-dimensional Ising and Potts models. Via this 
route the authors of also found expressions for one-time correlation functions of arbitrary order, 
valid however only for the dynamics after a quench from an uncorrelated initial state. 

In this paper we derive yet another version of the general solution for arbitrary correlation func- 
tions in the finite Glauber-Ising model subject to generic initial conditions. Compared to the existing 
results, our solution is in the most explicit form and cannot be simplified further for generic initial 
conditions. Like the authors of Refs. 0, @]> we start from the hierarchy of differential equations for 
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the correlation (or response) functions. However, instead of trying to solve the hierarchy directly as 
a whole, we judiciously split it into subsystems of inhomogeneous differential equations. These sub- 
systems are first solved separately; arbitrary correlation functions then follow in a recursive manner. 
Crucially, this recursion can be solved explicitly, giving in the end a rather compact solution of the 
full hierarchy This result is of interest in itself, but not directly useful for practical purposes since it 
involves summations over a number of terms growing with the system size TV. However, in the bulk of 
the paper we show that if the system is quenched from equilibrium at some initial temperature, then 
all sums involving an extensive number of terms can be performed. This gives the desired explicit 
expressions for arbitrary multispin two-time correlation and response functions, involving sums of 
only a finite number of terms. 

In Sec. [3 we start with a recap of the definition of the finite one-dimensional Glauber-Ising model, 
and give the hierarchies of differential equations governing the evolution of multispin one and two- 
time correlation (and response) functions. The general solution of these hierarchies, all of which 
have the same structure, is then derived in Sec. [H] We briefly discuss the procedure of taking the 
thermodynamic limit in our result, which is in fact straightforward, and then proceed with the rigorous 
analysis of the finite model. First we verify in Sec. IHII that the well-known equilibrium correlations 
at T > are retrieved correctly. From these, one-time multispin correlation functions after a quench 
from equilibrium at 7] > to arbitrary T > follow, as shown in Sec. II VI Here we make contact 
with the results given for the infinite chain given in [2], in ^| (infinite chain, 7] = oo and T = 0) 
and in (finite chain, 7] = oo). Using these one-time correlations as initial conditions we finally 
obtain the two-time multispin correlations (Sec. [Vj) and response functions (Sec. IVI|) after a quench; 
these are the core new results of our study. In Sec. I VIII we illustrate the procedure of extracting 
particular two-time functions from these general results. We briefly discuss the spin-spin correlation 
and response functions, already known from |lfl and then turn to domain- wall correlations 
and responses. There, as obvious applications, we give exact results and scaling expansions for the 
time-dependent equilibrium domain-wall autocorrelation at T > and for its two-time analogue 
for the out-of-equilibrium dynamics at T = 0. We also briefly discuss the connection between the 
Glauber-Ising chain and diffusion-limited reaction processes, indicating how e.g. multi-particle two- 
time correlations in the latter processes may be obtained from our results. We conclude in Sec. IVIIII 
with a summary and discussion of our work. Useful representations and identities for the functions 
appearing in our calculations are summarised in a number of appendices. 

Readers more interested in our results, rather than their derivation, could omit Sees. ITTHVII on 
a first reading. They need to refer only to the end of Sees. and IVII where we state explicitly 
the expressions for two and four-spin two-time correlation and response functions. The applications 
described in Sec. IVIII which are cross-referenced to the relevant sections in the derivations, should 
then make it clear how our general results are evaluated in practice. 



I. DYNAMICAL EQUATIONS 



The finite, one-dimensional Glauber-Ising model comprises a ring of N Ising spins a n € {— 1, +1} 
that evolve stochastically in time. We denote by p(cr,t) the probability to find the system in state 
<r = (eri, . . . (7^) at time t. Single spin-flip or heat-bath dynamics are described by the master equation 



d N 

— p(o-,t) = ^2 [w n (F n (r)p(F n cr,t) - w n (tr)p(tr,t)], (1) 

n— 1 

where w n (cr) is the rate for flipping the spin a n in state <r, and F n denotes the spin-flip operator 
F n <J = (<7i, ... — a n , . . . (Tat). Glauber's choice w n (cr) = i[l — ^a n {a n -i + oVi+i)] for the transition 
rates - with 7 = tanh(2J/T) - is such that the equilibrium solution of satisfies detailed balance 
and coincides with the Boltzmann distribution for the Ising Hamiltonian TL = — J^2 n a n a n +i at 
temperature T. The rates may equivalently be written as 

w ^ = 1 + c a!h W /t ' ( 2 ) 

where A n H(cr) = H(F n cr) — TL{cr) is the energy shift caused by flipping the n th spin. We remark 
that whenever site-indices exceed the range 1,2, ... N, as occurs for instance in the expression for the 
rates w±((t) and wn((t), we use iV-periodicity of the ring, that is <jn+i — <J\ and ctq = a N- 
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A. Correlation Functions 

We use the notation 

Cf } (t) = (an (*)••• a 4fc (t)) = 5><i • • ■ <r ijt p(<r, i) (3) 

for the (disconnected) fc-spin correlation function between sites 1 < i\, . . . ik < N at time t. The 
summation ^ ct obviously runs over the 2 N states of the system. From the master equation Q one 
readily obtains the evolution equation 

jj=1 

provided that the index- vector z = (ii, . has pairwise distinct components. The e v appearing 
in (@J are fc-dimensional unit-vectors along the ^-direction, that is ei = (1, 0, ... 0), e 2 = (0, 1, . . . 0), 
. . . efc = (0, 0, ... 1). Due to the index-shifts i ± e n , (@J induces a set of coupled differential equations 
over all fc-spin correlations. These equations are, however, not closed as for some % the index-shifts 
create a pair of equal components in i ± e n . In this case, the order fc of the correlation function is 
reduced by two, since <r„ = 1 for Ising spins; more precisely 

3 1 < M < v < fc : i„ = £7<*> ( t ) = C^. v) {t). (5) 

Here we have introduced the short-hand i\ (i^, i u ) = (i\, . . . i^-i, i^+i, ■ ■ ■ iv-\i iu+i, ■ • • ifc) to indicate 
that the components i^iv are to be removed from i. Since Q applies for any k > 1 and i with 
1 < ii, . . . ik < N pairwise distinct, the dynamics of (fc — 2)-spin correlations are again governed by 
J3J and in turn link to (fc — 4)-spin correlations, etc. For odd k this hierarchy of differential equations, 
linked by the hierarchical property (JSJ) of correlation functions, is closed at level k = 1 where the 
formation of index-pairs becomes impossible. The equations on even levels, on the other hand, are 
closed when an index-pair occurs on level fc = 2, giving C^(t) — °f Pi**-* *) = "l2a-P( (T ^) = 
Therefore the hierarchy of differential equations defined by Q , © over all even (odd) levels k' in the 
range < k' < k constitutes a closed set of equations for the evolution of fc-spin correlation functions 
with even (odd) fc, respectively. 

B. Two-Time Correlations 

We define multispin two-time correlation functions by 

C?f (t, O = K (t) • • • a ik (t)a h (*') • • • <7 i; (f )) = • • ■ a,, p(<r, i|<r', <') ^ • • • ^ p( < r', i')- (6) 

cr.er' 

Here p(cr,t\cr' ,t') denotes the conditional probability to find the system in state <r at time t > t' 
given it was in state a' at t'. This satisfies the master equation JIJ, 

d N 

— p(cT,t\a\t , ) = J2[MFn'T)p(Fn(T,t\a',t')-w n (a)p( ( T,t\(T , ,t')], (7) 

n— 1 

with the initial condition p(tr, t'\<r', t') = Sa-^r; here & a . a > = fl^Lj ^o-„,o-' is the Kronecker delta for 
states, which equals 1 for <j = cr' and vanishes otherwise. We deduce from J7J the evolution equations 
for two-time correlations 

|c^(m') = -kc% i \t, t >) + (c&J„,(*>o + ca2,j(*.o) • (8) 

77=1 

which again hold only for index vectors i with pairwise distinct components. It is clear from the 
definition that two-time correlations also exhibit hierarchical properties when index-pairs occur in 
i (or j). Therefore (SJ again induces a hierarchy of differential equations, linking (fc,/)-spin to 
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(k — 2, ^)-spin correlations etc. The hierarchies (JBJ) and (J3J have the same structure and differ only in 
the relevant initial values. For the two-time correlations these are 

c[f(t',t') = C^\t'). (9) 

The abbreviation i U j symbolises that the components of i and j are combined into a new vector 
(ii, . . . ik,ji, ■ ■ ■ ji)- The initial values © have to hold on each level k of the hierarchy for all i with 
pairwise distinct components. To avoid trivial reductions in the order I we restrict ourselves to vectors 
j with pairwise distinct components. 

C. Two-Time Response Functions 

In order to define response functions, we consider a perturbation STL — —Kb B{<t) to the Hamil- 
tonian, with a (generally time-dependent) field hs conjugate to the observable B(cr). The two-time 
response function RA,B(t,t') measures the response of the observable (A(t)) — A{cr)p{er : t) to a 
field impulse, where hs = h during a short time interval [t',t' + St'] and hs — at all other times. 
Transition rates in the presence of a perturbation are defined via J5J). Regardless of whether the 
system has been perturbed or not we may write 

(A(t)) = A((r) P (a,t\a',t' + St')p(a',t' + St'). (10) 

er.cr' 

From Q the conditional probability p(a,t\cr',t' + St') is unaffected by the perturbation since the 
field vanishes for t > t' + St' . To find the effect of the field on p(cr',t' + St'), one expands to linear 
order in St' to get 

N 

p{cr', t' + St') - p{a\ + St' Y, [w n (F n a')p(F n a', t') - w n (a')p(<r', t')} . (11) 

n=l 

Inserting this into (|TU)l and taking the difference between the perturbed and unperturbed cases gives 
the change 8(A{t)). Writing the expansion of the transition rates to linear order in the field as w n (a) = 
Wn(cr)\hB=o + hw' n {ar) then gives for the two-time response function RA^B(t,t') = TS(A(t)) / (hSf) 

N 

R A , B (t,t') =T^ A(a)p(a,t\a',t') £ [w' n {F n a')p{F n a' ,t') - W ' n {<r') p{<r' ,t')\ . (12) 

<r ,er ; n— 1 

Here an extra factor of T has been introduced into the definition of the response function, for later 
convenience. Since the derivation of equation Ijl2(l uses only the master equation Q), it holds for 
any choice of A(tr), B(a) and w n (cr). Comparing with © one sees that the ^-dependence arises in 
exactly the same way as for a correlation function between A(t) and some observable at time t'. For 
multispin response functions, where (A(t)) = (a^ (t) ■ ■ ■ c?i k it)) and STi. = —hj(t') aj 1 ■ ■ ■ aj x , one thus 
immediately has the dynamical equations 

|flft°(*.o - -*s8°(t.o +|E +fl2£j(t,o) ■ (13) 

ri=l 

These are to be interpreted in exactly the same manner as (@J) and ©■ To avoid a trivial reduction 
of the order / of the perturbation we require, in analogy to two-time correlations, pairwise distinct 
components in j. The initial conditions for l|13() at t = t' also follow from 1)12(1 : working out w' n (cr) 
from ((2J) gives 

i/=i 

Substituting into 1)12(1 . using p(cr, t'\cr', t') = Sa-^i and rearranging terms we then find 

fl—l U— 1 



<("•) = 



1 

2T 



-l^n+l) 



1 - 



c 



(k+l) 
iUj 



l_ r {k+l+2) 
2 iUj" 



(15) 
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In (I15JI . j v = (ji, . . . j v -i-,jv — lfjvtji/ + ■ ■ ■ ji) is the vector j with the additional components 

j v — 1 and j v + 1. The structure of (|15|l makes sense intuitively: an instantaneous response only 
occurs if the 'observed' sites ii, . . .ik overlap with the 'perturbed' sites ji, ■ ■ - ji- At level k = 0, the 
sum over [i in (|15Jl is empty; as is done conventionally, we define such empty sums as zero throughout. 
This makes sense: for k = 0, the observable is just a constant, A(a) = 1, and its value is unaffected 
by any perturbation, giving a vanishing response. Equations (@J) and JSJ are to be read in the same 
way for k = 0. 

II. GENERAL SOLUTION 

We have seen in Sec.[l]that the dynamical equations describing the evolution of multispin one- and 
two-time correlation and two-time response functions all have the same structure. Hence we obtain 
the various quantities from the same set of differential equations, e.g. I@J), by using the appropriate 
initial conditions. For this reason and in order to unify the analysis we use in this section the symbol 
pi k ) f or the dynamical quantity and for its initial conditions. For the sake of readability we omit 
indices i in the text except when they are needed for clarity. 

All dynamical functions F^ discussed above are by definition TV-periodic in the site indices it, . . .ik 
and invariant under permutations of the latter. But it turns out that these symmetries are not directly 
helpful for solving the hierarchy |@}. We therefore restrict the problem to functions F^ with ordered 
indices i £ N(k) where 

N(k) = {i e N k 1 1 < h < i 2 < ■ ■ ■ i k < N}. (16) 

The set N(k) is not closed under the index shifts i hi i - e, and i ^ i + e v occurring on the r.h.s. 
of as shifted components might exceed the range 1, ...JV or form pairs. Using JV-periodicity, 
permutational symmetry and the hierarchical property, however, each function F^ with a non- 
ordered index i ^ N(k) can be expressed in terms of another one, F^ k \ having an ordered index 
i' € JV(fc'). Hence we may formally rewrite (@J on all levels 2 < k < N as 

VieN(k): ±F?\t) = £ a^Ff{t) + £ &g,i^" 2) (i). (17) 

" j€N(k) j'eJVO-2) 

Note that (|17(l explicitly separates the terms on the r.h.s. of that link within level k of the 
hierarchy from those connecting to level k — 2 via the hierarchical property. Clearly the first term on 
the r.h.s. of Q, linking within level k for all i € JV(fc), produces diagonal entries in a^ k \ Also, each 
function with an index % — e v 6 JV(fc) and % + e n E N(k) on the r.h.s. of |@} yields corresponding 
off-diagonal entries in cS k \ But when a shifted index is not contained in N(k) several cases have to 
be distinguished. Here we only discuss the shifts i — e v as a similar reasoning applies to i + e n . If 
i — e v ^ N(k) and 2 < r\ < k the shift has created an index pair i^-x = i v — 1. Hence the hierarchical 
property applies and the corresponding term on the r.h.s. of (@J links to level k — 2. Such terms are 
accounted for in b^ k \ On the other hand, if i — e v ^ JV(fc) and r\ = \ then i\ = 1. JV-periodicity 
and permutational symmetry of the underlying dynamical function F^ allow us to replace the index 
i - ei = (0, i 2 , ■ ■ ■ ik) by (JV, i 2 , . . . i k ) and then by (i 2 , ...ik, JV), respectively. If (i 2 , . . . ik, JV) € JV(fc) 
there is a link within level k on the r.h.s. of (0J and we add a corresponding entry to in the 
appropriate place i = (1,42) ■ ■ - ik), j = (*2, ■ ■ - ik, JV). Otherwise ik — N is an index pair and this 
term produces an entry in . An explicit representation of as obtained by this construction is 
not necessary for our analysis; is given in (|A2|) for later reference. 

Rewriting (@J in the form Ijl7|l obviously does not change the basic problem of having to solve a 
hierarchy of differential equations. But l|17(l suggests a somewhat unconventional approach for finding 
F^: instead of simultaneously solving the full hierarchy induced by l|17l) we can focus on a given 
level k > 2 and consider (|17H as an inhomogeneous system of differential equations in F^ . Hence we 
may write 

t 

if) ( t ) = $f (t) + f At E - T ) b fj' F t 2) ^ ( 18 ) 

o jeN(k)j'eN(k-2) 
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where (t) is the solution of the homogeneous equation 

ViG A(fc): = £ ag^f (t), (19) 

satisfying the initial conditions Vi G A(fc) : ^^(O) = A^ k \ The second term in (|18|l . on the other 
hand, is a particular solution of (jl7(l obtained via the Green 's function of JTHJ). Thus we may 
first focus on solving the homogeneous equation (|19(l on the generic level k, obtaining the Green's 
function on the way. Then (|18|l constitutes a recursion formula that allows us to successively express 
any solution on level k in terms of solutions on lower levels k — 2, k — 4, etc. The recursion terminates 
either on level or 1 of the hierarchy where the solutions are known: according to the discussion in 
Sec.|I| for k = we actually have that F^(t) = AW is constant in time. On level k = 1, on the other 
hand, is homogeneous since no index pairs can be formed in a one-dimensional index vector. So 
FW(t) = is just the solution of (T^ on level k= 1. 

Note that the applicability of this approach essentially relies on the fact that links between different 
levels of the hierarchy are unidirectional, i.e. only from k to k — 2. This is a peculiarity of the one- 
dimensional Glauber-Ising model. 



A. Homogeneous Solution 

Solving the homogeneous equation l|19(l over ordered indices i € N(k) is a non-trivial problem. In 
order to make progress we apply the method of images, in this context already used by Glauber [Jj, to 
extend (|19|) over index vectors i € {1, . . . N} k . This procedure is purely for mathematical convenience 
and only possible because we have reformulated over ordered i G N(k) in the first place. As a 
generalisation of l|16[) we introduce the set of index vectors N„(k) that are ordered when permuted 
by 7T 

N n (k) = {i € N k 1 1 < Ml) < i <2) < . . . i <k) < N}. (20) 

Here, ir G S(k) denotes a permutation on {1, 2, . . . k} and S(k) the set of all such permutations. By 
(— l) n we refer to the sign of a permutation, that is (— l) p where p is the number of transpositions 
necessary to perform ir. Now we define the homogeneous solutions over {1,2,.. . N} k by filling the 
hypercube with permutationally antisymmetric images according to 

W e {l ) ...^:$f W = /(- 1 ) OT < ) W ^ W )(*) * 3*6 *(*):*€*,(*) (21) 

I otherwise 

This equation states that for ordered indices i € N(k) we have = & k K Furthermore, when i 
has pairwise distinct components any transposition in i just changes the sign of $( fe ). Therefore, and 
since $*- fe ^ = anyway for non-pairwise distinct indices, permuting the index i by ttq G S(k) changes 
the sign of <&( k ^ by (— l) 7r ° for any i G {1, . . . A} fc . Setting t = in (|21|l defines the initial conditions 
over the hypercube. As we explain below <J?^ satisfies 

k 

i* ■ A ^( fc ) m - _z~ m±lr fflsW ™ jl 



Vi G {1, . . . N} k : - ^ (t) = -fc + e„ (*) + Ce„ (*) J . ( 22 ) 

if the functions are moreover extended A^-(anti)periodically in each component of i for (even) 
odd k. We use the term A-antiperiodic in the sense that shifting any component of i by ±JV changes 
the sign of <£>( fc ). Note that in contrast to Q, which only holds for index vectors i with pairwise 
distinct components, l|22|l applies for all i G {1, . . . N} k . 

As a first step in proving that $( fc ) satisfies l|22|l we consider ordered indices i G N(k). According 
to (|21(l we have "l^ = for such i. Hence ()22l) should reproduce the homogeneous equation (|19|l . 
Equivalently, by expressing $( fc ' in terms of 'J^' on the r.h.s. of l|22|l we must recover the matrix 
a( fe ). In the first term on the r.h.s. of l|22|l . and in all terms with shifted indices i — e J( G N(k) and 
i + e v G N(k), we have <&( k ) — & k K So these terms directly reproduce the matrix entries discussed in 
the text below lfT7)l . Now we focus on the case where i — e v ^ N(k); as before, i + e r] N(k) is covered 
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FIG. 1: Illustration of the domain extension for homogeneous solutions in a system of size N = 6 and at 
level k — 2 of the hierarchy. The square in the center of the figure represents the region {1, . . . 6} 2 consisting 
of iV(2) (lower triangle) and its image (upper triangle). Associated with each circle in the lower triangle is a 
homogeneous solution $ = <I>. The upper triangle contains the permutationally antisymmetric image $ = — <f>. 
Along the dotted diagonal the index vector i = (11,12) is not pairwise distinct and hence $ = (full circles). 
The surrounding sections of squares represent the iV-antiperiodic extension of $ over Z 2 . Finally, the cross 
of arrows over a circle in position i points to the shifted index vectors i ± ei and i ± e.2 in correspondence 
with the terms on the r.h.s. of 12211 . Based on this simple example we consider the following special cases 
discussed in the text: (a) For i = (3,4) £ N{2) the shifted index vectors i + ei = (4,4) and i — e.2 — (3,3) 
are not in N(2) and must therefore contain an index pair (dots on the diagonal). The associated <E>'s are zero 
and do not contribute in J2Ub ( b ) For J = (2, 6) £ N(2) we have i + e 2 = (2, 7) N(2) (big circle). Using 
7V-antiperiodicity we replace $2,7 by —$2,1 which in turn satisfies $2,1 = —$1,2- Hence, this arrow effectively 
points to (1,2) £ N(2) and produces the correct term on the r.h.s. of 1221 1. (a',b') If 122H satisfies the cases 
(a) and (b) in N(2) it is immediately clear that it also holds for (a') and (b') in the image of N(2). (c) For 
i — (5, 5) we have e.g. i — ei = (4, 5) and i — e2 = (5, 4). Since $5,4 = —$4,5 these terms cancel in \2'll . Based 
on this illustration it is straightforward to prove cases like e.g. i = (6, 6) that are not discussed in detail in 
the text. 



by a similar reasoning. If i — e v <f N(k) and 2 < rj < k the shift has created an index pair = — 
According to (|21|) the corresponding is zero as required (case (a) in Fig.^). Alternatively, when 
i — e v ^ N(k) and 77 = 1 we must have i\ = 1 and therefore i — e.\ = (0, £2, . . . ik) ^ {1, • ■ ■ N} k . Using 
A?"-(anti)periodicity we replace this index by (N, 12, ■ ■ ■ ik) € {1, ■ ■ ■ N} k and modify the sign of <J>W 
by (— l) fe_1 . Next, applying a cyclic permutation tt c turns this index into (12, ■ ■ .ik,N) and changes 
the sign of $^ fc ' according to lf2*T|l once more, by (— l)*" . But (— l) 7r ° = (— l) fc_1 in fact cancels the 
sign-change caused by 7V-(anti)periodicity. So, all in all, if (£2, ■ ■ - ik, N) € N(k) we recover the correct 
matrix element of between i = (1,^2, • • - ik): 3 = (*2, ■ ■ • ik, N) (case (b) in Fig. 0. Otherwise 
ik = N and hence $' fe ' = as required. 

Next consider indices i £ N v (k) for some given ir € S(k) that are ordered when permuted by 
tt. Here we use the fact that is antisymmetric under permutations for all i € {1, • • .N} k . So, 
regardless of whether a shifted index on the r.h.s. of (122[l is also contained in N„(k) or not, we permute 
it by 7T and change the sign of <&( k ) by (— l) 71 ". In cases where a shifted index lies outside {!,... N} k 
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we first use A-(anti)periodicity, then permute by tt and finally apply 7r c . As discussed above, this 
again just gives a sign change (— iy . Due to linearity of (|22|) the overall factor of (— l) w drops out. 
Therefore tt maps each sector i £ iV ff (fc) of ill' 21) onto % e N(k). Since 122(1 holds on i £ JV(fe) this 
implies that 1(22(1 is also satisfied over i G N w (k) for any 7r S S'(fc) (cases (a',b') in Fig. GJ. 

We complete the proof that <&( fe ) satisfies ((22(1 over i € {1, . . . A} fe by showing that it holds for 
non-pairwise distinct indices, too. In this case there must be at least one index pair, say = i„, in 
i. Thus the l.h.s. of 1(22(1 . the first term on the r.h.s. and each term with a shifted index 77 5^ ^i,v 
contain this pair and are identically zero according to 1(21(1 . We also know that swapping the /I th and 
v th component of i — or i + turns it into i — e u or i + e„, respectively, since = i v . But it 
takes an odd number of transpositions to swap index-components, consequently the associated $( fc )'s 
have opposite signs and cancel each other in 1(221) (case (c) in Fig. GJ . When index shifts exceed the 
hypercube {1, . . . N} k we use iV-(anti)periodicity and the same argument applies. 

Summarising so far, by appropriately extending the definition of d?^ 5 ) beyond the physically relevant 
range i S N(k) we have managed to rewrite the homogeneous equation ((19(1 in the much simpler form 
1(22(1 . A comparison of ((22(1 with the original hierarchy J3J shows that antisymmetry has naturally 
eliminated all links between different levels k of the hierarchy. Now we can solve ((22(1 over the space of 
./V-(anti)periodic functions. For this purpose we introduce the discrete Fourier transform f q = J~{f n } 
as 

N 

VneZ: f n = ^Y.fi e+mq - ( 24 ) 

If one restricts 1 < n < N in ((24(1 . then ((23(1 , 124(1 form a transform pair when Q is any equidistant 
TV-partitioning of the interval [0, 2tt]. But in order to retrieve iV-periodic or A-antiperiodic functions 
from 1(24(1 the particular choice Q — Q or Q = Q e has to be made, respectively, where 

= {f + f n\n = 0,...N-l} = {§,%,...&^} (25) 

Q Q = { %n |n = 0,...A-l} = {0,f,...^^}. (26) 

So on even levels k of the hierarchy we set Q = Q c as we are dealing with A-antiperiodic functions 
while Q = Q on odd levels corresponds to TV-periodic functions. In either case the transforms (|25jl . 
124(1 satisfy the usual shifting and convolution properties of Fourier transforms. Hence the fc-fold 
transform of 1(22(1 follows immediately as 

Vq E Q k : 1 (t) = (t) + 7 £ cos ' ( 27 ) 

Here q = (31, • • • Qfe) are the Fourier variables corresponding to i = . . . Note that we use the 
subscript to distinguish $^ from its transform $g fe ' . In Fourier space the evolution equations 1(27(1 for 
<l>(' c ) are decoupled and can therefore be integrated easily. Then, by inverting the Fourier transform 
and using the initial condition $^ fc ' ) (0) = we obtain 

*S fc) (*) = E ne-^^^^lf, (28) 
h,—jk »7=1 



where the sum runs over j S {1, . . . A} fc . An expression for I n (x) is given in I(B1(I in Appendix IbI 
where we also discuss the properties of I n (x) that are relevant for our analysis. Finally, by expressing 
the initial conditions in terms of A^> via 1(21(1 and after some rearranging, we obtain the solution 
$( fc ) of the homogeneous equation ((19(1 on level 1 < k < N 

VieN(k): *<*>(t) = E Gf}{t)Af (29) 



jeN(k) 



with the Green's function 



Vi>0: Gg(t) = E ("irfle-* 4,-^(7*)- (30) 

7re5(fc) r,= l 
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The Green's function as such is defined over the index range i,j £ N(k). Expression H3()|) . however, 
carries the antisymmetry of the images and is therefore well defined over i,j e Z fe . For later conve- 
nience we note that (|30|) is in fact permuationally antisymmetric in i and j. To see the antisymmetry 
in ji, consider the sum over tt e S(k). If, instead of j, we substitute a version permuted by ttq, we 
may use tt o ttq G S(k) as the summation variable since the permutations form a group. Hence the 
replacement of j by (i^m; ■ ■ -j-K {k)) just changes the sign of the Green's function by (— l) 7r °, prov- 
ing permutational antisymmetry. Antisymmetry over i can be shown similarly and is an immediate 
consequence of the images. 



B. The Recursion Formula 



In this section we use the recursion formula (|18|1 to construct the general solution of the hierarchy 
l|17|) . It is convenient to express the homogeneous solution $w contained in (|18|) in terms of the 
Green's function G^> via l(2*9l . This gives 



= E G?}{t)Af+ UtY, E Gf}{t-r)bfl,F^\r). (31) 

jeN(k) I jeN(k)j'GN(k-2) 



Now we show that a particular sequence of simplifications, essentially based on the similarity of 
different levels of the hierarchy, can be applied when 131JI is iterated. Our proof is inductive, hence 
we start with proposing a structure for the result. Working out the first few iterations of (|3*T|) suggests 
that L + 1 < Lf J iterations should give the following expression 

^i fe) w = E E (-irn^^(»)C^.) (t)+ < (i)< (32) 

1=0 w6 p(; i fc) A=l 

Here V(l,k) denotes the set of permutations that correspond to choosing < I < |_§J ordered pairs 
from the numbers {1, 2, ... k} and keeping the remaining k — 21 numbers in ascending order, i.e. 



V(l,k) = {tt € S(k) | tt(1) < tt(2), tt(3) < tt(4), . . . ir(2l - 1) < 7r(2Z), 
tt(1) < tt(3) < . . . tt(21 - 1), n(2l + 1) < tt(2Z + 2) < . . . Tr(k)}. 



(33) 



An expression for H n (t) is given in IjGljl . We also discuss relevant features of this function in Appendix 
ICl The appearing in (|32|l are the homogeneous solutions l|29|l . Finally, denotes the remainder 
term of the recursion that links to level k — 2L — 2, and is of the form 

/drl E (-irn^»A)-M«-i,(2(*-T)) (34) 

{7r£V(L,k) X=l 

x y y G ( ( k - 2L) . , .(t-r)b (k - 2L) F^- 2L - 2 \r) 

j£N(k-2L) j'eN(k-2L-2) 

The first step in our proof that l|32l) is the product of iterating (|31() is to consider L — 0. The two 
sums in Q32J) collapse to a single term, that is I = and "P(0, fc) = {Id}. This obviously also applies to 
the first sum in l|34|l . Defining the value of empty products in (|32|l and (|34(l as one - a convention we 
use throughout then reduces (|32|l to l|31|) for L = 0. Now assume that (|32|) is true for some L > 0. 
Consequently, for i|32|) to apply also for L — > L + 1, the remainder has to satisfy 

7r£P(L + l,fc) A=l 

Of course we need L + 1 < [|J - 1 since (|31|l only holds on levels fc > 2. We show below that the 
remainder H34JI satisfies H35JI which proves that (|32J) is the result of iterating (|31(l . 
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First we evaluate the matrix product G^ k 2L ' • 2L ' in (|34|) : we show in Appendix El that this 
can be expressed in terms of Green's functions G( fc ~ 2i ~ 2 ) and the functions H n {t). Substituting (|A1|) 
into EH 

gives 



= l dT { E E (-i)^- 1 

■rr£V(L,k) l< fJ .<v<k-2L 

L 



-°-H 



(36) 



E 



G 



(fc-2L-2) 



(V(2Z,+1),— V(A,))\(V(2L+ /4 ),V(2t+v)),i 



j'€JV(fc-2i-2) 



Then we express 

F (k-2L-2) in J2gj 

via the recursion formula lj31|l and use the general identity for 

our Green's function 



Vr,r'>0: £ Gg (r) Gf}, (r') = G$(r + r'). 
jeJV(fe) 



This gives 



(37) 



ir£V{L.k) 



l<H<v<k-2L 



■H, 



A=l 



E 

jeN{k-2L-2) 



/•»' s 



(» w (2L + l),---« T (fc))\(«^(2i + M ),V(2L + ^)).J V 7 J 



(fc-2L-2) 



E 



jGN(k-2L-2) j'£N(k-2L-l) 



G 



(fc-2L-2) /, _ /» 

(V(21 + l).---i7r(fc))\(»-7r(2i+ M ).»- n -(2I, + i/)),i^ ' J). 



,-.Ak-2L-2) „(k-2L-4), / 



V) 



(38) 



In l|38|l r-dependencies via the Green's functions have dropped out. 

Next, we rewrite the permuations in l|38[) in terms of those from the set V(L + l,fc). Consider 
the case L > 1 first. According to l|33|) the first 2L values of tt e V(L,k) correspond to L ordered 
pairs chosen from {1, . . . k} whereas 7r(2i + 1), . . . n(k) are the remaining numbers in ascending order. 
Now the /i, v sum chooses another ordered pair ir(2L + fj,) < tt(2L + v) since /i < v and tt is ordered 
in that range. In contrast to the L pairs selected by tt the new pair does not obey the mutual 
ordering between pairs. As both sums run over all possible choices we get each set of L + 1 pairs 
L + 1 times, with ir(2L + fi) in all mutual orderings from tt(2L + /i) < 7f(l) < . . .tt(2L — 1) to 
7r(l) < . . .tt(2L — 1) < tt(2L + /i). But this ordering is only relevant for the second line in (|38ll . 
Hence we may restrict the choice of pairs to those contained in V(L + 1, k) and confine the sum over 
the L + 1 ways of selecting a "special" pair to the second line in (|38[) - which thereby turns into the 
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derivative of the product of L + 1 terms 



l_7reP(L+l,fc) V A=l 

G (fc-2(L+i)) 



.-(«_!) (2(*-r)) 



> .ft) A. 



(fe-2(L+l)) 



(39) 



jeN(k-2(L+l)) 



/•'••• £ 



jeAT(fc-2(L+l)) j' S JV(fc-2(L+l)-2) 



G 



(fe-2(L+i)) fi-r'U^" 



M,(fe-2(i+l)) „(k-2(i+l)-2)/ , 



3 



if) 



Now we convince ourselves that the sign of the permutation in (|39|l matches the sign- factor in (|38[1 : 
in H38[) we start from the ordered vector (1, 2, . . . k). Applying it S V(L, k) permutes this vector as 
described, with the sign of the permutation being (— l) 71 ". The sum over \x, v then chooses components 
2L + /i, 2L + v as pair number L + l. Now it takes v — /i — 1 transpositions to move those components 
next to each other and an even number of further ones - thus irrelevant - to move this pair around. 
Altogether, (— l) 7r (— l) ly_M_1 in l|38|) is the sign of the corresponding permutation 7r 6 V(L + 1, £;) in 
(El- 

In the L = case, equivalence of (13811 . (|39|l is obvious: since 1^(0, k) = {Id} we can remove this 
sum and all 7r's from (|38|l : the empty product over A is one by definition. So l|38(l is just the explicit 
formulation of the sum over it £ V{1, k) in (33$ . 

Finally an integration by parts in t turns Ij39(l into (|35|l . using H n (0) = 0. This completes the 
proof that (|3"2|) is the result of iterating (|3"T]) . 

Having proved (|32ll . we can now set i = |_§ J — 1- The remainder term in (|32ll then links to level 
or 1 of the hierarchy for k even or odd, respectively. But, as discussed below (|19[) . we know the 
functions = and F^ = $W appearing in the remainder. In fact it can be simplified in full 
analogy to the above calculation for L < [f J — 1 and this yields the general, explicit solution of the 
hierarchy (|17fl as 



E E 



( 1) ^i,r(2A)-M2A 

2=0 TreV(l,k) A=l 



,J2t) <S> { , k ~ 2l) . , 



(t). 



(40) 



The homogeneous solutions <&( fe ' appearing in l|4(JI) are given by (|29[l for 1 < k < N. We have 
incorporated the link to level in (|40l) . which occurs for even fc, by defining — A^°>. The definition 
of the sets V(l,k) may be found in (|33|l and expressions for the functions I n (x) — contained in the 
$W — and H n (t) are given in l|Bljl and IjCljl . respectively. According to H4()|l all functions I n (x) 
and H n (t) appearing in are associated with (even) odd levels for (even) odd k and consequently 
the appropriate sets Q = Q e or Q = Q as given by Eqs. I|25I26|) must be substituted in their 
representations. Explicit versions of (|40[1 for the first few levels read, bearing in mind that the 
product over A equals one for I = 0, 



*f>(t) 



*< W (t) 



Ff (t) = <f>( 2) (t)+i/ 42 _ 4l (2t)$(°) 



ff>(t) 



+ ff <a _ <1 (2t)*W(i) - ff <8 _ <1 (M)*^(t)+fli,- <a (2t)*^(i) 



F i2 _ il (2t)$g )i4) (t) -^ i ,_ il (2t)$w <4) (t) + 5^(2*)*^ 
H i3 _ i2 (2i) $g, i4) (t) - H^ l2 {2t) *g A) (t) + H u _. l3 (2t) d>g >i2) 
f 44 _ 43 (2t) - fffe-H (2t) F i4 _ i2 (2t) $W 



(2) 



Hi 2 -i 1 (2t) ^ 



i (2t)$ (0) 
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We remind the reader that the functions only have direct physical meaning if the indices i are 
chosen from N(k). Extended over all indices i G {1, . . . N} k they are permutationally antisymmetric; 
this follows from l|31[> and the corresponding property of the Green's function. Therefore the at 
any time t can be used directly as antisymmetrized initial conditions A^ k > in (|28J) : this property will 
be useful below. 

Equation (|40|l is evidently the simplest possible representation of F^ for generic initial conditions 
as encoded in the $( fc )_ While for the generic case the complexity of F^ k ' grows rapidly with the 
level fc, we will see below that in physically interesting cases significant simplifications occur. 



C. The Thermodynamic Limit 

The entire analysis presented in Sec. [H] can be repeated in full analogy for an infinite spin chain. 
There, one takes the thermodynamic limit in which thereby becomes an infinite hierarchy of 
differential equations. The individual levels k may again be rewritten in the matrix form i|17[) over 
ordered indices. For the infinite chain, however, one chooses — oo < i\ < i% < . . . if. < +00 to obtain 
complete equations. The infinite matrices and b^ k ' have a slightly simpler structure since the 
shifts i — ei and i + ej, can never produce an index pair. In other words iV-periodicity of the finite 
ring is lost. Nevertheless (II Til may be treated as an inhomogeneous system of differential equations 
which can be solved according to the steps presented in Sec. Ill Al and Sec. Ill Bl Rather than the 
discrete iV-(anti)periodic Fourier transforms (|23l) . (|24|l the standard Fourier series expansion is useful 
for solving the homogeneous equations. The structure of the solutions 140|) of the hierarchy and the 
homogeneous solutions (|29|) remains completely unaffected. But instead of the functions I n (x), H n (t) 
their N — * 00 limits l n (x), H ra (<) discussed in Appendix IbI ICl appear in Ij29(l . (|40|l . The sum over the 
initial conditions in l|29|) obviously runs over — 00 < ji < . . . % < 00 for an infinite chain. 

Alternatively, the rigorous solutions (|29|l . (|40() for the finite ring also produce the same N — > 00 
limit. This is evident for H40[l as it comprises a finite number of functions, the limit of each of which 
exists. In (|29|l . however, one has to keep in mind that the functions I n (x) are iV-(anti)periodic. 
Therefore all summation indices j that are close to i modulo N contribute to the sum in (|29l) : the 
site N for instance is a neighbour of site 1. So a meaningful N — > 00 limit can only be taken if we 
first shift the index range of the finite ring. An appropriate replacement for N(k) is for instance 
— [ r 1 j < i\ < ...ik < LtJ- Shifting the index range on the ring changes the domains of the 
solutions (|29p. 14U|I but leaves them unaffected otherwise. Convergence for N — > 00 of the sum in 
l|29(l over the symmetric index range can then be proved |l6j . for any fixed index vector i. 

In the following we continue to analyse the finite ring of spins. But according to the above discussion 
the thermodynamic limit may be taken at any stage. It just amounts to replacing the functions I n (x), 
H n (t) by their N — > 00 limits I„(x), H n (i) and extending summation ranges. From now on we omit 
the limits for summations that are taken over the ring, as for instance in l|28|l . In a finite system 
these summations run over 1, . . . N while in the thermodynamic limit they are to be evaluated over 
—00 . . . + CO. 



III. EQUILIBRIUM CORRELATIONS 

As discussed in the beginning of Sec. |n] the functions F^ (t) — (t) describe the evolution of 
correlations if we use the initial conditions = C^ k \0). Here we show that our solution l|40|) 
reproduces the well known equilibrium correlation functions for T > 0. 

Equilibrium correlations arc obtained from the dynamical solutions Q40(l by taking the limit t — > 00. 
Let us first focus on the homogeneous solutions $( fc ). It is convenient to consider the representation 
J2EJ instead of J35J|; recall that for i S N(k) we have $ (fe) = $ (fe) . When taking the modulus of (25)1 
and using the triangular inequality we may drop the initial conditions < |C^ fc ^(0)| < 1. This 

gives the first inequality in 

*f(*)|< E n e_t 1 J *n-^ C-y*) I = II e ~* 53 1 I < e-**/--*. (41) 

jlf-jkV^ 1 '7= 1 J-n 

As indicated, the fc-dimensional sum in j|41|l factorizes once the initial condition terms are removed. 
Then, using the bound (|B10I) on J2 n gives the second inequality in (4"T)) for any system size 
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N. Here l/r cq = 1 — 7 is the smallest eigenvalue of the master operator corresponding to 0). 
Since T > implies 7 = tanh(2 J/T) < 1 the equilibration time r eq is finite. Therefore at any level 
k > 1 of the hierarchy the homogeneous solutions vanish in equilibrium. At level k = 0, however, we 
have $(°)(t) = A(°) = = 1 at all times. 

Next consider the i — * 00 limit of the functions H n (t). Based on the representation (|C2I) it is 
straightforward to show that for T > and arbitrary system size N we are left with the sum 

ff„, cq = lim H n {t) = ^Y] sm(n<?) 7Smg . (42) 
t— 00 iV ^— ' 1 — 7COSO 

By combining the results ijSJl and (|42|l we now obtain the equilibrium correlations from (|40|l . 
Any correlation function C' ( - 2 ' c+1 - ) containing an odd number of spins is expressed in terms of 
$W , $( 3 ) , . . . $( 2fe + 1 ) . But all these homo geneous solutions vanish in equilibrium and hence so do 
the correlations of odd order, as expected. For equilibrium correlations of even order C^ 2k ^ , on the 
other hand, the only nonzero terms in 140JI are those having a factor $' ' = 1. This gives 

He 

c t k i = £ (-irn^M^-M--^ ( 43 ) 

TreV(k) \=1 

where we have introduced the notation V(k) = V(k, 2k) for the usual set of ordered pairings. Since 
the functions H nteq in l(43() are associated with even levels we need to set Q = Q e when substituting 
1(421) . For ordered indices i <G N(2k) the result 1(43(1 describes multispin equilibrium correlations in the 
finite Glauber-Ising model at T > 0. Over non-ordered indices i € {1,. . .N} 2k , however, (|43|l carries 
the permutational antisymmetry of 1(40(1 : this property will be useful below. 

The link between 1(43(1 and the result obtained from a transfer-matrix calculation is established by 
realizing that for 1 < n < N — 1 we have for (|42|l the identity 

1 ^ . . . 7sin(7 [tanh( J/T)} n + [taah(J/T)] N ~ n 

— > sin tio — ■ — = - — — v ' n . (44) 

N ±1 ^l- 7 cosg l + \ta,nh(J/T)} N 



96Q 



The l.h.s. of (|44|l is identically zero for n = 0, N and A^-antiperiodic, however. To prove (|44|) one 
notes that the sum may be rewritten in the form of the inverse Fourier transform (|24l) . So we may 
equivalently prove the converse statement Ij23(l for q e Q c - The latter just involves a geometric 
summation and is straightforward. We finally remark that when the representation 1)44(1 is used for 
H n , C q, the permutational antisymmetry of (|43|) is lost. Only for ordered index vectors i S N(2k) are 
all combinations n — i^^X} ~ V(2A-i) m ESt i n the range where l|44() applies. 



IV. CORRELATIONS AFTER A QUENCH 

Next we derive dynamical correlation functions C^ k \t) from our general solution (|40|l . We consider 
the example where the system is quenched at t = from an equilibrium state at 7; < 1, corresponding 
to T; > 0, to some arbitrary temperature 7. As in Sec. IIIII we may identify F^ k '(t) = C^ k '(t) if we 
use equilibrium correlations as the initial conditions A^ k \ 

On odd levels of the hierarchy we have A^ 2k+1 ' = since equilibrium correlations between any odd 
number of spins vanish. Consequently $( 2fc + 1 ) = according to lj2"§j) and via also C( 2fc+1 ) (t) = 
at all times. This conclusion is trivial because the master equation and the equilibrium initial 
state are invariant under spin inversion er — > cr. On even levels 2k > 2 of the hierarchy we express 
A (2k ^ via g2J and (gSJ) as 

A m = £ (-D* ft £ sin [ q ( M2A) - M2A _ 1}) ] T ^£_. (45) 

As pointed out below 1(43(1 . this representation for the equilibrium correlation functions is permuta- 
tionally antisymmetric and consequently zero for non-pairwise distinct indices i. This is why we have 
been able to write A instead of A in ((45(1 . Knowing A is convenient as it allows us to derive the 
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homogeneous solutions from 1(28(1 instead of the more complex expression ()29(l . In fact, since (|45J) is 
factorized over index pairs, the sum in 1(28(1 also factorizes into two-dimensional sums 

1 • m 7ising 



E e 2t ^(2A)-m(7*)^ ( 2A-i)-™(7i) E sin [« ( m _ n )l i 



7i cos q 



(46) 

Using the properties (EJl, (|B8|) of the functions I n (x), the sum ^ m n in (|4*6*|) can easily be evaluated. 
Here it is essential that the functions I n (x) in 1(46(1 are associated with the even levels 2k and that in 
the argument of sin[g (m — n)] we have q £ Q e too, as (|B7|) would not apply otherwise. We obtain 

*? h Ht)= E i-^UK^-^-^ ( 4? ) 

where H' n (t) is the result of the summation over m, n and may be written as 

Kit) = 1 E sin(ng) 7iSlng e-^-v-^). (48) 

We remark that this straightforward derivation of the homogeneous solutions heavily rests upon 
the permutational antisymmetry of our expression 1)43(1 for equilibrium correlations. When Glauber 
derived the dynamics of two-spin correlations he expressed the equilibrium values via ((44(1 . He 
introduced images at this stage of the calculation to make ((44(1 permutationally antisymmetric. 

Having worked out the homogeneous solutions ((47(1 that correspond to our problem, we now obtain 
the correlations after the quench from 140(1 

cf k \t) = E E (-ir E II K^-^m n ^U+^vb-^^v-,^ 2 *)- 

2=0 7r£P(Z,2fc) fteV(k-l) A=l A' = l 

(49) 

For the term with I = k to correctly give the contribution from & k ^ — = 1 one has to define 
■p(O) = {Id} here. To simplify 1(49(1 . we now focus on the combination of pairings that occur for some 
fixed I: each tt £ V(l, 2k) draws / ordered pairs from {1, . . . 2k}. These define the indices of the factors 
H n (2t) in ((49|) . Out of the remaining, ordered 2k — 21 numbers n(2l + 1), . . . 7r(2fc), the permutation 
7T £ P(k — I) selects another k — I ordered pairs, which give the indices of the factors H' n (2t). If we 
denote by 717 £ S(2k) the overall permutation on {1, . . . 2k}, then all 717 satisfy 

7T Z (1) < 7H(2), 7T Z (3) < 71,(4), . . . Wl(2k - 1) < 7T Z (2fc) 

(50) 

tt;(1) < 7r ; (3) < ...7^(21- 1), tt ; (2Z + 1) < n (2l + 3) < . . . 7n(2* - 1). 

By ordering the first I pairs against the subsequent k — I ones, each 717 £ S(2k) can be mapped onto 
an ordered pairing tt £ V(k). This leaves the sign of the permutation unaffected as it takes an even 
number of transpositions to reorder pairs, i.e. (— 1) T = (—1)"' = (— l) w (—1)*. Conversely, we have 
that for each tt £ V(k) there are ( ; ) semi-ordered pairings 717 of the form 1(50(1 . These correspond to 
the different ways of distributing k — I primes among the functions H n (2t) in 1(49(1 when we choose 
the index pairs according to ft. Together with the summation over I we can thus factorize ((4911 for 
each ft, giving 

cf k \t)= E i-^UK^-^-^ ( 51 ) 

where H^(t) = H n (t) + H' n (t). Since the functions H n {t) in ((49(1 are associated with even levels 2k 
we need to set Q = Q e in the representation l(C2p . Together with ((48(1 this yields 



= J- y S in(n g ) J 7i Smg e-^-T""^ + 7Smg 1 - e^ 1 -""" k (52) 



-t(l — 7 cos g) 



TV ^ \ 1 — 71 cos (7 1 — 7 cos (? 

At t = our solution for dynamical correlations J5TJ, ((52(1 obviously reproduces the equilibrium 
values 1(42(1 . ((43(1 for the initial temperature 71. Note that for a quench from T; = 00 (71 = 0), i.e. 
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a random initial configuration, l|52|l reduces to H^(t) — H n (t). It is equally clear that we recover 
the equilibrium values at 7 in the limit t — > 00. The spectrum of relaxation times for correlations 
is encoded in the exponents in (|52() and depends on the final temperature 7 after the quench and 
the system size via Q c as given in l)25[l. In the thermodynamic limit N — > 00 we can replace ^2 q 
by ^ J dq in 152|) and, for consistency, then use the symbol H'^(i). Further properties of H^(t) are 
discussed in Appendix [0 An equivalent formulation of the results l|51|) , l|52(l for the infinite spin 
chain was derived in 0]. 

It is worthwhile to note that the general solution (|40|l preserves the structure of ordered pairings. 
In the case at hand, for instance, we used the initial conditions <|43|) and arrived at the result (|51|l . 
Both expressions have the same structure and (|51(l again carries the permutational antisymmetry 
of l)4()|l. Therefore the dynamical correlations after an arbitrary sequence of temperature shifts are 
still of the form (|51l) and could be written down immediately as long as the system is initially in an 
equilibrium state. Examples of two and four-spin correlation functions as obtained from (|51J) read 

cf\t) = K^t) 

Cf ] (t) = (2i) Hl_ l3 (2t) - //;' , (2t) <_, 2 (2i) + Hl_ h (2t) H^_ %2 (2t). 

The number of terms grows as the number of permutations contained in V(k), i.e., (2k — 1)!! = 
1 ■ 3 ■ • ■ (2k — 1), giving 15 terms for C^ & \ 105 terms for and so on. 



V. TWO-TIME CORRELATIONS AFTER A QUENCH 



As discussed above the evolution equations (jHJ of two-time correlations C^ k ' 1 ^ (t, t') may be written 
in the form l|17|) . Hence these functions again follow from the general solution H40|) if we use the 
equal time initial conditions 10 at t = t'; obviously time arguments t in (|40|l have to be replaced by 
A = t — t 1 in this setup. We consider the same quench as in Sec. lIVI because this allows us to express 
the initial conditions (© in terms of (|51|l for even k + 1. On all levels k where k + 1 is odd, that is all 
even levels k for odd I and vice versa, we have C {k ^(t',t') = C {kJrl \t') = 0. Hence C {k ^(t,t') = at 
all times from Ij29(l . I|4Q(I : this is again just a consequence of the spin inversion symmetry er — * — cr. 
So from now on we assume that k + I is even. 

In analogy to Sec. II VI an expression A for the equal time initial conditions 10 that is permutationally 
antisymmetric in i would be desirable. Without loss of generality we require that the vector j - 
only appearing as a parameter — is ordered, j S N(l). Now the following problems have to be 
handled. First, if i € N(k) is also ordered, our expression for A has to produce the correlation 
function between the k + I sites i U j . But the combined vector is not necessarily ordered and hence 
the associated correlation cannot be expressed directly via l|51jl. Second, there could still be pairs of 
equal components between i and j. While the correlation over i U j then reduces its order via the 
hierarchical property JSJ, the expression (|51|l is identically zero. A systematic construction of A based 
on @ and l|51(l is possible but rather cumbersome. Let us therefore immediately state the result for 
A and prove its correctness instead: 



\{k,l)(.i 



( II W (iUi) w(2A _ 1) ,(iUi)„ (2A) (*') 



(53) 



where 



I] sgn(j A - Q 

A=l 



11 sgnljA ~Q sgn(jA - is) 



K-tJV) 



(a,b) = (juJv) 
for (a,b) = (i e , j v ) (- 54 ) 
(a,b) = (i e , is) 



The sgn(n)'s in (|54|l are the standard sign functions, with sgn(n) = — sgn(— n) = 1 for n > and 
sgn(0) = 0. To illustrate the notation of l|53|l and l|54|l we consider as a simple example the 6-spin 
two-time correlation with k = / = 3. In this case the sum in l|53|) runs over the 15 ordered pairings 
contained in V(3). Among these focus on the particular pairing ir that gives (1, 3, 2, 4, 5, 6). Here the 



1G 



indices (a,b) of the 3 factors U in iftEfy are ((tU (iU j)s) = (ii,i3), ((iUj)a, (*Uj)4) = («2,ji) and 
((*Uj)5, («Uj) 6 ) = {j2,h)- According to 0j2J we substitute n A s S n 0'A -«i)sgn(j A -^3)^"-^ (2t') for 
the first factor Wj 1 ,i 3 as both indices are drawn from i. The second factor I4i 2 j 1 comprises components 
from i as well as j and must therefore be expressed as (—l) 1 ~ 1 Si 2t j 1 + Y[\ s & n (j>> — ^)H" i _ i2 (2t l ). 
The third factor, with both indices from j, is simply Uj 2 j 3 = H" a _^{2t'). The remaining terms in 
the sum over tt G V(3) in (|53l) are worked out similarly. 

A few further comments on (|54|l arc in order. First, the case (a, b) — (j„,i e ) cannot occur in 1(54(1 
because all pairs in tt are by definition ordered. Second, the W's carry the superscript j since the 
products in (|54|l run over the full vector j, so that U depends on all components of j except in the first 
case of <|54|) . Finally, the factor (— 1)" _1 in 1(54(1 depends on the actual number v of the component 
j v drawn from j . 

Having clarified the interpretation of ((53(1 and (|54|l we now prove that these equations are indeed the 
initial conditions A corresponding to ©■ Let us first show that our expressions are antisymmetric 
under permutations in i € {l,...iV} fe . It is sufficient to prove antisymmetry under transposition 
of any two adjacent components of i, say 1^,1^+1, since all permutations may be decomposed into 
such transpositions. Now, given i^i^i, we split the pairings V in 1(53(1 into those V 1 containing 
the pair (i^z^+i) and the remaining pairings V" where i M and are distributed over separate 
pairs. For each tt G V' there is a factor Ui i +1 in ((53(1 which, according to (|54|l . is of the form 
n A sgn(j A - in)sgn(j x ~ i |t+1 )fl-^ +i _ i(i (2f). With H'l n (t) = -H%(t) this immediately shows that 
Ui = —Ui i +1 ; therefore each term in (|53|l with tt G V is antisymmetric under «-> 

For tt S T 5 " we have to apply a different reasoning. Suppose, for instance, that a particular pairing 
tt G T 5 " contains the pairs (i e , i^) and As the pairs are ordered we have e < fi, 6 < fi + 1 

and the mutual ordering e < 5. But 5 ^ [i implies S < /i and hence the conditions e < fi + 1, 
S < /i and e < d also hold. So there is another pairing n' G V" that contains (i £ ,i M+ i) and (isji^) 
but is identical to 7r otherwise. Since tt and 7r' have opposite signs and because the transposition 
V V+i maps the pairings corresponding to tt and 7r' onto each other, the sum of the two associated 
terms in l|53|) is antisymmetric under <-> i^+i. This argument is readily extended to all tt G V" 
by considering the remaining cases; these are pairs of the form [i E , i^), (i^^x) where e < \x and 
x G {v+2, ■ --ik,h, ■ ■■3l] or (ifj,,x), with x, y G {v+2, • --ji}- 

Now we turn to ordered index vectors i G N(k). Then (|53|l must produce the correlation function 
over the sites i U j according to the initial condition 0. Let us for the moment focus on those 
i G N(k),j G N(l) for which the combined vector iUj has pairwise distinct components. This allows 
us to ignore the Kronecker delta in 1)54(1 . Thus U 3 a b contains a factor J\ x sgn(j\ — i e ) whenever the 
index a is drawn from i, i.e. when a = i e , and similarly for b. Since all components of i occur exactly 
once in the product over A in (|53"|) we can therefore rewrite (|53"|) as 

k I — 

Apart from the products over the sgn's this corresponds to our result 1(51(1 for correlation functions. 
But ((51(1 is antisymmetric under index permutation. Hence, if tt G S{k + I) is the permutation that 
orders iUj, i.e. iUj G N^(k + l), then the pairing sum in ((55(1 must give (— l) w times the correlation 
over xUji. Thus the products over the sgn's must coincide with (— l) w in order to cancel the sign of w. 
We can convince ourselves that this is indeed the case. Keep [i fixed and evaluate the product over 
v. For each v with j v < we have sgn(j„ — i M ) = —1 while z M < gives sgn(j„ — i M ) = +1 and is 
irrelevant for the product. Altogether this gives (— l) rlfl where j„ < < j n Since both i G N(k) 
and j G iV(Z) are ordered we have n\ < n.2 < • ■ ■ ^fc- Now if we start ordering i U j beginning with 
ifc, it takes transpositions to turn i U j into (ii, . . . ife— 1, Ji, ■ ■ ■ j« fc , ik,jn k +i, ■ ■ - ji), a further rik-i 
transpositions to turn this into [ix, . . .i fe _ 2 ,ji, . . ■ jn k -i,h-i,jn h ^ x +i, ■ ■ ■ jn k , ik, jn k +i, ■ ■ ■ ji) and so 
on. Therefore the products over the sgn's give (—1) to the power of the number of transpositions 
necessary to order i U j, which is equal to (— l) 7r . 

Finally we lift the restriction that i U j must have pairwise distinct components; i G N(k) and 
j G N(l) is of course still required. Assume that there are p pairs i £l — js 1 ,i S2 — js 2 , ■■■is p = js p 
between i and j. We write k' = k - p, I' = I - p, i' = i \ (i El , . . .i £p ) and j' = j \ (js 1 , ■ ■ -js p )- 
The combined vector t U j then contains no further component pairs. According to the hierarchical 
property J^J the equations fiufy and must then reduce to (55) with k,l,i,j replaced by k',l', 
respectively. The first step towards showing that this is true consists in splitting the pairings in 
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again. Now, however, we denote by V' the pairings that contain all the pairs (i £q ,js q ) with q — 1,. . .p 
and by V" the remaining pairings. Let us focus on tt £ V" first. For each such tt at least one pair, 
say (i £q ,js q ), is not formed. Instead i £q may be paired up with (i u ,i £q ), (i e , i v ) or (i e ,j w ) where 
w 7^ S q . In either case there is a factor U = in the product in 1(53(1 since Ui^,i F , U% E ,i„ and 
all contain a factor sgn(jj — z e ) =0 according to 1(54(1 : for Ui c j w the Kronecker delta is also zero 
since j £ N(l) is ordered and hence j w ^ js — i £q for w ^ 8 q . So each term in the sum over tt in l(53|) 
is zero for tt £ V" . Next we turn to the pairings tt £ V' ■ Each such tt may be constructed as follows. 
Starting from iU j = (ii, . . .i £l . . .i e2 . . .i £p . ..i k ,ji ■ - .jst • •• j<s 2 -• -js p ■ ■ ■ jz), we first group all pairs 
(i £q ,js q ) with q — l,...p into the lowest components of the combined vector. It takes k — e\ + Si — 1 
transpositions to move js 1 next to i £l and an even (irrelevant) number of further transpositions to 
get (i El ,j Sl ,i 1 ,...i £2 ...i £3 ...i £p ... i k ,j 1 ...jg 2 ... js 3 . . . js p ■ ■ ■ ji). Moving js 2 next to i £2 then only 
takes k — £2 + 82 — 2 transpositions because component jg 1 is no longer contained between i £2 and j$ 2 
in the latter vector. Clearly, regrouping i U j into the form i rg = (i £l , , i £2 ,js 2 ,..- i £p , js p )Ui Uj 
then requires a permutation 7r rg of sign 

p 

= JJ(-l)*- e «+ tf «-«. (56) 
9=1 

Now by arranging the top k' + I' components % U j of i rg according to tt £ V((k' + l')/2) and 
ordering the pairs (i £q ,js q ) against ((*' U j')*(2A-i), (*' u J )*(2A))j any pairing tt e 7"((fc + Z)/2) 
of i U j may be obtained. Also, since rearranging pairs takes an even number of transpositions, 
(— 1)* = (— l) 7r ' E (— 1)^. Using this representation for the pairings tt £ V' we may now rewrite 15 311 . 
First we substitute the factorisation for (— 1)^ just given. Second, since all pairings tt contain the 
pairs {i £q ,js q ) with q = l,...p each product in i(5*5|l has the factors Mi e ,j Sq which just give (— 1) A « _1 

according to 1)54(1 . The remaining factors are over the pairings % of i U j . Since % U j has pairwise 
distinct components we may apply a similar reasoning as in l(55ll and thus obtain 

^SV) = (-i)^n(-i)**- 1 n ri^-v) 

fi^ei,...e p 



k' + l' 



X 



E ("I)* II ^W)^,-^')^^^')- (57) 



From 1(56(1 we see that equation (|57|l reduces to (|55|l and therefore implements the hierarchical property 
(JSJ for pairs between i and j if the identity 



9 =1 



(U#ei,... 



(58) 



holds. The product over /x is easily evaluated when using that i £ N{k) and j £ N(l) are ordered 
and i £q = ja . Hence sgn(j$ q — i^) = — 1 for fx = s q + 1, . . . k. Since /1 = £ q +i, . . . £ p are omitted, the 
product over /j, gives (— l) k ^ £ i-p+i _ So the factors in the product over q become (— l)~ 1-p , yielding 
the overall product (— 1)~p(p +1 ) = 1 since p(p + 1) is even for any integer p. This completes our proof 
that equations 1)53(1 and 1(54(1 are the initial conditions A corresponding to (0. 

Based on the expression (153(1 for A it is straightforward to write down the corresponding homoge- 
neous solutions $ via 1(28(1 . Every component of i in 1(53(1 — which becomes a summation variable in 
(|2*%|l — occurs in exactly one of the factors of the product over A in ijftlfy , paired up with (depending 
on the particular pairing tt) either a different component of i or with a component of the fixed vector 
j. Once the fc-dimensional summation over i from ((28(1 is exchanged with the summation over tt from 
l(53() . it thus factorizes into one- and two-dimensional sums. This procedure yields 

k + l 

*$ o (t,o= e (-^nw-.w),,,,^'), (59) 
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where V denotes the case-dependent sums over U 



K b (t,t')= { (-l)"- 1 e- A J i ._ J>( (7A)+^ ijV (t^) for (a, 6) = (i e , j„) (60) 



' J ,,--W \^ _i_ £?' 

fi, is (t,t') (a,b) = (i e ,i s ) 

and £, T are, respectively, the one and two-dimensional sums 

dim(j) 

S L»W) = EI1 S gn(jA-m)e- A / le _ m ( 7 A)^_ m (2i') (61) 
m A=l 

dim(j) 

rLisQJ) = E II sgn(jA-m)sgn0 A -n)e- 2A / Je _ m ( 7 A)/ 45 _„( 7 AX'_ m (2t'). (62) 

m,n A— 1 

Equations (|59J) and (|60|1 are to be read in exactly the same manner as Ij53(l and Ij54(l . The sums over 
m and n in i|61|) and (|62|l each run over TV different values, so that these expressions are not practical 
for large system sizes. We show in Appendix [Dj however, that both sums can be reduced to ones 
involving only a finite number of terms, so that the limit TV — > 00 can be taken without problems. In 
<|61|) and l|62l) we have replaced the dimension I of ji, being the upper limit of the products, by the 
generic expression dim(j). This is for later convenience. Note that according to i(55|) and l(6"UI) the 
homogeneous solutions for arbitrary order two-time multispin correlation functions can be expressed 
purely in terms of the functions /, H" , £ and T . Examples of low order homogeneous solutions read 

$ ( °' 2) (M') = K-nW) (63) 

= e- A / 4 _,( 7 A)+^(M') (64) 

S>?f(t,t>) = Fl M {t,t')Hl_ h {2t>) (65) 
- [ + c- A I ix _ h ( 7 A) + £l n (t, f )] [ ~ e- A I^ 32 ( 7 A) + B^ (t, t')] 
+ [ - e- A ( 7 A) + Bl h (t, f )] [ + e- A I^ 3l ( 7 A) + B?^ (t, t')] 

Substituting the result (|59|l into (|40|l finally delivers analytic expressions for the desired two- 
time correlation functions. We have worked out C^(t, t'), C<: 2 > 2 ) (t,t'), (t,t'), C^ 3 ' 1 ) (*,*') and 
C^ 3 ' 3 ^(t, t') and it appears that a factorisation similar to the one which simplified 149JI to (|51Jl is 
generally possible. We do not have a proof for this claim, however, and therefore take (|40|l together 
with H59f) as our final result, which expresses an arbitrary multispin two-time correlation function as a 
sum over a finite number of terms. Examples of low order multispin two-time correlations, expressed 
in terms of homogeneous solutions via (|40|l , read 

c\Y\t,t') = ®\f\t,t') (66) 
cg 2) (M') = <&g 2) (t,^)+^- 1 (2A)$5 2) (i,^) (67) 

In Sec. IVIll below we present a discussion the correlation functions C^ 1,1 ) (t,t') and C^ 2:2 \t,t'). In 
particular, we demonstrate how the expressions l|66|l and l|67|l simplify when using the results of 
Appendix IdI to rewrite the sums iffiTJl and (jH2l- 



VI. TWO-TIME RESPONSES AFTER A QUENCH 



Quite analogously to the case of two-time correlations discussed above, two-time multispin response 
functions j?( fe '') (£, t') follow from the general solution (|40|) with time arguments replaced by A = t — t' . 
The equal time initial conditions R^ k ' l '(t\ t') are given in terms of correlation functions in (|15fl . We 
consider again a quench from equilibrium as introduced in Sec. IIVI Consequently, on all levels k 
with odd k + I the initial conditions (|15|l vanish and so do the R^ k ' l \t,t') for all t > t' . Thus we 
immediately focus on k + I even. Without loss of generality we restrict the vector j, labelling the 
sites to which the perturbation couples, to ordered j S TV(Z). 
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Equation 1(15(1 comprises a correlation function over the sites i U j", where j" is the index vector 
j with the additional components j v — 1 and j v + 1 as defined below l|15|l . Since j is fixed we can 
immediately use the hierarchical property (0 of correlation functions to remove component pairs 
from j v in cases where j v £ N(l + 2). More precisely this corresponds to replacing j" in ((15(1 by the 
"squashed" index vector 



(jl, ■ 


■■ju-l 




,ju 


Ju - 


I" 1 Ju+1, ■ 


■ ■Jl) 


(?!.■ 


■ ■ ju-2 




Ju 


Ju - 


Hi Ju+1,- 


■ ■Jl) 


{jl, ■ 


■■ju-l 


Ju-l 


Ju 




7 ju+2 , ■ 


■ ■Jl) 


0'i, • 


■ ■ ju-2 




Ju 




, ju+2 , ■ 


■ ■Jl) 



for 



ju-l <jv-l,jv + K ju+1 
ju-l = ju -l,jv + l< ju+1 
ju-l < ju -l,j„ + l= ju+1 
ju-l = ju -l,j„ + l= ju+1 



(68) 



We denote the dimension of j v ' s by l v ' s € {I — 2, 1, 1 + 2}. Having made this modification in 1(15(1 . which 
guarantees that ji"' 55 6 N(l"' s ) is ordered, we may now use (|5rt(l to write down the initial condition 
A corresponding to l|15|l : recall that 1(53(1 is the permutationally antisymmetric version of ©. This 

gives for the permutationally antisymmetrized initial condition for (t, t') 



k I 



- EE*~ 



1 i 

„2 



1 ~ y j e ( _i ) 7r n w (<Ltf)«(»-i),(<ui)«(») c* 7 ) 



2 e ( _ n 



A=l 



A=l 



i ),r(2A-l),(*Uj , '- s )„ (2A) 



(0 



(69) 



The fact that the Kronecker deltas in (|69|l contain components of i does not affect permutational 
antisymmetry: the only occurrence of the summation variable /i in (|69(l is in 8i j v . But finite sums are 
independent of the order in which they are taken, hence permutational antisymmetry is unaffected. 

We now note that for each /i, v the Kronecker delta in ((69(1 . enforcing z M = j v , allows us to simplify 
the expressions inside the square brackets. As discussed in the text below Eqs. 1(531 154(1 . the products 
over A are zero for all pairings that do not contain the pair (i^Jv). So for each /i, v it is sufficient to 
restrict the sums over the ordered pairings in 1(69(1 to those that do contain the pair (i^Ju)- For the 
pairings drawn from iUj, the corresponding factor U with index (a, b) — (i^Ju) is (— l) 1 ' -1 according 
to But for the pairings drawn from iUj"' s this factor is — (— l) u 1 since j v is component number 

+ 1 or v — 1 of j u ' s as can be seen from 1(68(1 . Thus we may rewrite ((69(1 as 



/i= 1 V— 1 
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7_ 
2 



E II W '( J iUj)x(2A-l):(iUj) x(2A) 



(0 



A=l 



E ^ 1 ) ,r II W '( J iUj>'. 8 ) ir(2X _ 1) ,(iUi''. s )„ (2A , 



(0 



A=l 



(70) 



Here 7-" emphasises that we only sum over the pairings that contain {i^J v ) for each given fi, v. 
Since we have pulled the factors ±(— \) v ~ x — contained in each product in (f?0"|) — out of the square 
brackets we have to replace U by 1 for (a, b) — [i^Ju), and indicated this by writing W instead of U. 

From the permutationally antisymmetric initial condition ((70(1 that corresponds to ((15(1 one easily 
obtains the homogeneous solutions for two-time multispin response functions via ((28(1 . We use the 
fact that according to our construction of ((70(1 the entire expression inside the square brackets is 
independent of i M for each fixed fi. So the summation over this index is trivial. The sum over the 
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remaining indices i \ (i M ) follows in full analogy to the derivation of 1)59(1 and consequently we obtain 



fl—l u—l 



( 7 A) 



2 ^ 



)^(2J.-l)>( i Uj"' B ) 7r(2 A) 



(71) 



Here V" = 1 for indices (a, 6) = (i^,j v ) and V" = V as given by ((60(1 otherwise. As in the case of 
two-time correlations we see that the homogeneous solutions for arbitrary two-time multispin response 
functions can be expressed entirely in terms of /, H" ,£ and T according to 1(71(1 . At level k = 0, the 
sum over fj, in (|7T)l is empty and so we set $(°'') = A^ '^ — as discussed at the end of Sec. II CI In 
Appendix El we explicitly demonstrate how the first nontrivial homogeneous solutions (|7S|l and f74l 
are obtained from the general result ((71(1 : 



$(0.2) 
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(72) 
(73) 
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^2-ii+i 



-fii-ii+l + £ 



- t »2— .72 — 1 i °; 



O'l-l Jl) 
«2 J l-1 

(ii-i,ji) 

iijji— 1 

0'2J2+1) 



r 1 c Oi J2 i 



J2 — 1 



«2j'2 + l 
C 0'2,J2 + 1) 



(74) 



To save space we have omitted the arguments of the functions / and £ which are obviously 7A and 
(t,t'), respectively. The expression 1(74(1 for $( 2,2 ) only holds if j is of the form j = + 1); for 
ji + 1 < j2 the result iffTjl gives a more complicated form. 

As usual, actual two-time response functions may be expressed in terms of homogeneous solutions 
1(71(1 via 1(40(1 . It is not clear to us whether a factorization of the combined expression 1(40(1 . 1(71(1 is 
possible. So we state our results explicitly in terms of homogeneous solutions, as for instance 



RffW) = <f>i 2 - 2) (t,i')+^2-u(2A)$ ( "' 2) (t,^) 



(75) 
(76) 



The responses R^ 1 ' 1 '^,^) and R( 2 ' 2 } {t, t') are discussed further in the subsequent section. In partic- 
ular, we show that the identities for £ given in Appendix IdI yield a simplification of the expression 

Ci- 



vil. SIMPLE EXAMPLES OF TWO-TIME FUNCTIONS 



In this section we illustrate the procedure of extracting explicit expressions for multispin two-time 
correlation and response functions from the general solutions derived in Sees. IVl and I VII For some of 
the examples we also give a discussion of their implication for the physics of the Glauber-Ising chain. 

The simplest possible two-time correlation functions comprise only one spin at each time. By 
combining ((64(1 , ((66(1 and rewriting the sum £ using l(D6(l we obtain 

C\y\t,t') = = e-^frAJ + fl^AX). (77) 

Expressions for the functions I n (x), H^(ti,t2) are given in l|Bl(i . ((Clip . Since the two-time correla- 
tions 1(77(1 are the solutions of JSJ) at level k = 1, the set Q = Q a must be used in (|B1|) . The result 
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1)77(1 applies for the finite Id Glauber-Ising model quenched from an equilibrium state at some 7] > 
to an arbitrary temperature T > at time t — 0. According to the discussion of the functions I and 
H in Appendices [B] and [CJ taking the thermodynamic limit N — > oo and quenching from a random 
initial configuration T; = oo just amounts to replacing Ij-i and H'-_ i in 1(77(1 with and Hj-j, 
respectively. Writing out (A, 2i') explicitly using I(C11(I we thus have 

C\f\t, = e- A Ii_,-(7A) + sgn(j - i) H^A, 2f ) + 7 / drc^ I^t) Hi (A - r, 2*')- (78) 

Jo 

If we additionally quench the system to T = then ((77|) together with the zero temperature formula 
l(C23() for E n (ti,t 2 ) simply yields 



hi-j(t + +jf drli- j(t + f - r) [I + Ii] (r)| 



Cg' 1) (t ) t / )=e- (t+ *'M^ ■,(/ + /')+/ <lrl, ;(/-/'-D [T„ + I,.](r) }. (7<i) 

Here we have introduced the short hand [ • ] (x) to indicate that all functions contained in the square 
bracket have the same argument x. The functions l n (x) in 1(78(1 . ((79(1 are the modified Bessel functions 
(jB2j). 

Similarly we obtain the two-time spin response functions to a local magnetic field hj in the non- 
equilibrium state after the quench - scaled by T according to our definition in Sec. II CI - by putting 
together (Eft, G5J 



(80) 



I„(a;) and H' 7 [{t) are given in I(B1(I (with Q = Q Q ) and ((52(1 . Taking the thermodynamic limit AT — > 00 
and letting Ti = 00 turns ({HQ} into 



4\' 1) (^0=e- A Ii- J -(7A) 



y(l + H 2 (20) 



(81) 



For the zero temperature quench 181|) may be simplified further by using l(C17(l . giving 



RffM = l - e-( t+t ') I 4 _,(A) [Io + 2I X + l 2 ] (2t>). (82) 

The exact solutions 1(78(1 ■ 1(81(1 agree, as they should, with their temporal Laplace (with respect to 
A, t') and spatial Fourier (with respect to n = i — j) transforms as derived in [l(]j. While the latter 
are surprisingly simple for the functions ((78(1 and ((81(1 . a corresponding transformation of higher 
order functions is virtually impossible. The correlations ((78(1 and a modified version of the responses 
((81(1 . slightly differing in the way the transition rates in the presence of a field are defined, were also 
studied in |lll | ; the long-time scaling form of the correlation functions can equivalently be found using 
a picture of annihilating random walkers • 

Both publications [ToL Hl| considered the non-equilibrium fluctuation-dissipation relations asso- 
ciated with the auto-correlation and response functions. For any pair of correlation and response 
functions one can define a fluctuation-dissipation theorem (FDT) "violation" factor X via 

R(t,t') = X(t,t')^C(t,t') (83) 

so that X(t,t') — 1 corresponds to the usual equilibrium FDT (recall that our response function 
contains an extra factor of T; in the conventional definition, l(83|) would have X/T on the r.h.s. instead 
of X). Under appropriate conditions, T/X can then be interpreted as an effective temperature T e fj ; 
this interpretation is well-established for mean-field spin glass models but its extension away from 
mean-field remains open and has prompted much recent work |l8l Il9l |20| . 

In we have given a detailed discussion of the FDT violations that result for non-local correlation 
and response functions, extending the work of [Tfl.lTll | for the local quantities. Here we therefore only 
illustrate one aspect of our results, by considering the local quantities obtained for i = j from (|78^l . 
I(81|) for a quench to nonzero temperature. When all times t, t' and r cq = 1/(1 — 7) are large, one 
expects results that depend only on the time ratios t/t' and i'/T eq ; this is indeed what we find, in 
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agreement with the results given in . A convenient way of representing the results is a parametric 
fluctuation-dissipation (FD) plot of the integrated response function against the correlation function; 
the former is defined as 

X (t,t') = J drR(t,T) , (84) 

and gives the response to a field switched on at t' and held constant after that. One notes from 
<|83|) and R(t,t') = —(d/dt')x(t,t'), that the negative slope of such a p lot is guaranteed to give 
X(t,t') only if t' is used as the curve parameter while t is held fixed [21|. This gives the FD plots 
in Fig. fright); on the left we show the FD plot reported in which was produced with the 
opposite convention of keeping t' fixed and varying t along the curves. We see that the plot on the 
right produces a physically far more plausible representation of the crossover from the aging regime 
t' <C T cq to equilibrium t' ^> r cq ; in particular, one reads off that the FDT violation factor X(t,t') is 
always < 1 and its value in the limit C — > evolves smoothly from 1/2 to 1 as the ratio t/r cq grows 
from to oo. 




C(t,t') C(t,t') 

FIG. 2: FD plots of the local spin-spin correlation and response functions, for fixed t' (left) and fixed t 
(right). The values of t' /r cq are t' /r cq = 10~ 3 , 1CT 2 , 10~\ 1CT 1/2 , 10°, 10 1 from bottom to top on the left; the 
same values are used for t/r oq on the right. See text for discussion. 

Among the higher order correlation functions, one of primary interest is the domain wall correlation. 
The domain wall indicator (or defect variable) rii = \ (1 — <Xj erj+i) G {0, 1} measures the presence 
(rii = 1) of a domain wall between sites i and i + 1. It is a well known fact [ESQ that Glauber 
dynamics of the (unmagnetised) spin system <x corresponds to the diffusion limited reaction process 
A + A ?= in terms of n — (m, . . .n^); here we interpret rii as the particle occupation number 
of site i. The relevant rates follow immediately from (J2J): the transitions ... T J, J, ... ^ ... tTJ, •■■ 
for spins, where | and | symbolise 0{ = +1 and Ui = —1, respectively, occur with rate 1/2 and 
correspond to diffusion ... 10 ... ^ ... 01 ... of particles. The process . . . |J,| . . . ^ . . . flf . . ., on the 
other hand, which maps onto ... 11 ... ^ . . . 00 . . ., gives the rate (1 + "f)/2 for particle annihilation 
and the (1 — 7)/2 for pair creation. For zero temperature dynamics, where 7=1, this obviously 
reduces to diffusion limited pair annihilation (DLPA). In the context of the Glauber-Ising model, 
understanding the dynamics of domain walls or defects is crucial for e.g. explaining the coarsening 
process after a quench. Here our general solutions for multispin two-time correlation and response 
functions given in the preceding sections are the key for a comprehensive analysis. Due to the mapping 
to diffusion reaction processes, our results simulatenously describe multi-particle two-time correlation 
and response functions in the latter systems. 

As a simple example we now consider the two-time connected correlation function between two 
domain walls, or equivalently particles, Ci-j(t,t') = 4[(n,i(t) rij(t')) — (n,i(t)) (rij{t'))]. In terms of 
spins this gives 



Gi-jtyJ) = (tr i (t)cr i+1 (t)cr i (t / )o- i+ i(0) ~ M<) <?i+i(t)) {a 3 {t') a ]+1 {t')) . (85) 
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Now consider the coarsening dynamics after quenching the system from an equilibrium state at T\ > 
to some temperature T > 0. The two-spin correlations in l|85|) follow immediately from (|5ip. i.e. 
(o"i(i) (Tj+i (i)) = H'{(2t). The non-trivial term in (|85() is the two-time four-spin correlation function. 
Combining the results ((63(1 . (|(>f>|) and 1)67(1 we may write any such correlation as 

Cff{t,t') = [i/ 42 _ 4l (2A)+^^ 2 (i,0]^" 2 -, 1 (2t') (86) 
- [ + c- A 4 ( 7 A) + el dl (t, t')] [ - e- A I i2 . h ( 7 A) + £l j2 (t, 0] 
+ [ - c- A 4 _ J2 ( 7 A) + ^ (t, f )] [ + e- A I i2 . h ( 7 A) + f/; J . i (t, f )] . 

To obtain the two-time correlation in the first term of (|85|l we set i — (i, i+ 1) and j = (j, j + 1). 
For these index vectors and bearing in mind that ff" = —H'_!_ n , we obtain from l|D3|) and (|D8|) the 
following representations for the sums £, T in ((86(1 

= ^'_ i (A ) 2i')+e- A / i _ J _ 1 ( 7 A)^'(2t') 

f/; ja (t,0 = i?;_ i+1 (A ) 20-e- A / i _ J -( 7 A)Hl'(20 

^(t.f) = i?;_ i _ 1 (A,2t')+e- A 7^( 7 A)^(2i') 

El h (t,t') = ^"_ 4 (A,2i')-c- A /^ +1 ( 7 A)iJr(2t') 

= ^(2A,2f) 

- e" A [/^ ( 7 A) flJLi+i (2*0 - (7A) Hti (2*')] 

- e" A [/ i _ J -_ 1 ( 7 A)^ J .(2t / )-J < _ J -( 7 A)^L J -_ 1 (2t')] 

+ e- 2A / 4 _,( 7 A) 7 i _ J ( 7 A) JTf(2f ) - e-^J^-xC-yA) I t - 3 +i( 7 A) iJf(2t') 

Substituting these expressions into 1(86(1 . using the identity (|CT0)l to rewrite H n (2A) + #"(2A, 2i') = 
H"(2t) and cancelling terms gives 

cQ 2) (t,t') = H'{(2t)H'{(2t') (87) 
+ [e- A /^( 7 A) + fl/'_ i (A,2f) ][e~ A /, w ( 7 A) - ffj'_ 4 (A, 2f) ] 
- [ e- A J^-xCtA) -fl/'_ m (A,2f)] [ e- A ( 7 A) + fl/'_,_ 1 (A,2f )]. 

Note that the first term on the r.h.s. of (|57|). #j'(2t) H'{{2t'), cancels with (cr 4 (i)cr l+ i(^)) (<7j (*')) 
when we substitute into (|85|l . The correlations 1(871) are those in a finite ring of spins; since 1(871) was 
obtained as the solutions of © on level k — 2 the set Q = Q e has to be used in the representation 
JEJ for I n (x). The functions H^(ti,t 2 ) are given m {C9jl . Now we focus on the quench from a 
random initial state T; = oo and take the thermodynamic limit -/V — > oo, thus replacing I n (x) and 
H^(tx,t2) in l(87|l by I„(x) and H„(ti, t2)j respectively. Substitution into 1)85(1 then yields the two-time 
connected defect correlation function 

C n (t,t') = [e~ A I„( 7 A) + H n (A,2f) ][e" A I„( 7 A) -H n (A,2t') ] (88) 

- [ e- A I„+i( 7 A) - H„ +1 (A, 2tf) ] [ e- A I„_i( 7 A) + H n _!(A, 2t') ] . 

For the T — quench ((88(1 may be simplified further using l(C17(l 

C n (t,t') = e-( H -^p n _ 1 -I n+ i](t + t / )Hn(t-t / ,20 (89) 
+ e- 2 *I„(; - i')Pn-i + 2I n + Wi](* + - e- 2(t+t,) [(I„-i + I„)(I„ + l n+1 )](t + t'). 

For small n the functions H n (ti, £2) in ((89(1 may be expressed purely in terms of modified Bessel 
functions via the recursion ltC17|) . For large n, on the other hand, it is more convenient to use 

the representation l(C4(l . with I n {x) replaced by l n (x) and 7 = 1. 

We have used the zero temperature result ((89|l in ^1 to study the non-equilibrium fluctuation- 
dissipation relations for defect observables n^. Beyond this, we are not aware of any expressions 
equivalent to e.g. (|8*%|) in the literature. As explained above, our results are of interest not only for 
the coarsening dynamics of the Glauber-Ising spin-chain, but also directly give the corresponding 
particle-particle correlation and response functions for the associated diffusion reaction process. 



24 



From H88J) we can, for instance, derive an exact expression for the equilibrium domain- wall autocor- 
relation. This function corresponds to (four times) the particle autocorrelation in a diffusion limited 
reaction process A + A =± at equilibrium. It follows from (|88|l by setting n = and taking the limit 
t' — > oo with A = t — t' fixed. This gives, using the representation (|C4|I for H„(A, 2t'), 



{/>oo 
e- A Ii( 7 A)-J/ dre^[I -I 2 ]( 7 T) 
1 J A 

Expanding asymptotically for A, r e 



oo at fixed A/r Gq gives the scaling forms 
2 

1< A < r cq : C eq (A) 
C cq (A) 



A 



1 < r cq < A : 



' c q „-2A/t„ 
ttA 2 



(90) 

(91) 
(92) 



Plots of H90(l for various r cq are shown in Fig. In the regime A <C r oq the particles may be 
considered as an ensemble of independent random walkers. Hence the decline of C eq (A) corresponds 
to the return probability of a random walker. In the opposite regime, however, the particles are 
extremely likely to have been annihilated and replaced by new ones via the process A + A t— 0. Since 
these new particles are uncorrelated with the original ones, the connected correlation vanishes. 

For comparison we also consider the non-equilibrium case of diffusion limited pair-annihilation 
A + A — > 0, for which the correlation function is given by the zero temperature formula l|89[) . Setting 
n = again yields an exact expression for the particle autocorrelation 



C (t, t') = 2 e" 2 * I (A) [I + h](t + - e- 2 (*+*') [I + h] 2 (t + 1') 



(93) 



In H93() the time t' plays a role analogous to r cq in H9()(l . So we expand for A, t' — > oo with A/t' fixed. 
This gives the scaling forms 



1 < A < t' 
l«i'«A 



C (t,t') 
C Q (t,t r ) 



2? 



irA 2 



(94) 
(95) 



Plots of (|93|l are also shown in Fig. [21 For A <C t' we have a situation that is exactly analogous 
to the equilibrium case, with T eq replaced by 2-7rf'. In the opposite regime A 3> t', however, the 
correlation function decays much more slowly, now due to the annihilation of particles that are 
very far apart rather than, as in equilibrium, the creation of new particle pairs. We discuss in |T^ 
that the scaling form l|95() is nontrivial and, for instance, gives rise to unanticipated noncquilibrium 
fluctuation-dissipation relations in the Glauber- Ising model. 

The counterpart of the two-time connected domain wall correlation function H85JI is the two-time 
domain wall response function to local bond perturbations hjj+i defined by 



R^ 3 {t,t')=T 



S(a t (t) a l+1 (t)) 



Shj, j+1 (t') 



1=0 



In full analogy to the correlations we obtain R n (t,t') from l|74() . 17611 and l|D3(l . with i 
j = (j.j + 1). For N — > oo and a quench from a random initial state T; = oo we find 



(96) 

(i, i+ 1) and 



R n (t, t') 



e" A I„ 


[-(' 


e- A I„-i 


K 1 


e" A I„+i 


H 1 


e- A I„ 


[-(■ 



T, 



1-1 (-e- A I„ 



f+e- A I„_! 



T, 



H„ + i 
H n _i 
H„ 



^(+e- A I 

A 



—A 



)-!( 



-1-71+ 1 



H„_ 2 j 

H n _iJ 
H„+i) 



^ (-e" A I„ +2 + H n+2 ) 



(97) 



All I„ in (|97|l have argument I„(7A) and all H„ = H„(A, 2t'). For the T = quench we may use the 
recursion (|C17J| for H„(£i,t2) to simplify this equation drastically. It turns out that the result may 
be written as 



d 



Rn(t, = ^ {e" 2t I n (t - [I„_i + 2I n + I n+1 ] (t + t')} 



(98) 
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FIG. 3: Connected domain wall autocorrelation functions Co(t,t') = 4[(ni(t)m(t')) — {ni(t)){rii(t')}] in 
equilibrium (left), i.e. t' oo, and for the coarsening dynamics after a zero temperature quench at t = from 
a random initial state (right). The different curves correspond to decreasing temperature (left) and increasing 
system "age" %' (right) from top to bottom at small A = t — t'. The plots are obtained by numerical evaluation 
of Eqs. (19UI93I for r cq , t' = 10 1 , 10 2 , . . . , I0 6 and are discussed in the text. 



The functions C n ({55j) and R n il5Yjl are rather special: for other index values, e.g. + 2), or 
in particular for higher orders k, I the expressions become rather more complicated and are most 
efficiently evaluated using symbolic software. As above, simplifications occur for N — > oo and for a 
quench from 7] = oo to T = 0. Any two-time multispin correlation and response function C^ k ' l \t, t') 
and R^ l \t,t' ) with k and I even can then be expressed purely in terms of modified Bcssel functions 
l n {x), due to ljHT7)l . (|D3)l and JUHJ). It is clear from (|Hl7j) . of course, that the number of I n (a;)'s in 
the result grows with the index range covered by 



VIII. CONCLUSIONS 



We have presented a new approach for solving the full hierarchy of differential equations for mulit- 
spin correlation and response functions in the finite one-dimensional Glauber-Ising model. Our result 
is the most explicit representation for these functions. The known results for equilibrium and dy- 
namical correlation functions as well as two-time spin-spin correlation and response functions after 
a quench are easily recovered from our solution. Beyond that, however, we have derived closed ex- 
pressions for arbitrary mutispin two-time correlation and response functions for the dynamics after a 
quench. We found that any such quantity can be expressed purely in terms of the four functions I, 
H" , £ and T ', regardless of the number of spins involved. 

Our results for the two-time correlation and response functions formed the basis for the study Ijjj , 
where we gave a comprehensive analysis of the nonequilibrium fluctuation-dissipation relations in 
the Glauber-Ising model at zero temperature. A number of other applications of our results can be 
envisaged. For example, based on the general expression for two-spin two-time correlations, we are 
currently studying the existence of dynamical hereogeneities in the Glauber-Ising model |22| . 

Finally, we have emphasised the link between the Glauber-Ising chain and diffusion-reaction pro- 
cesses. In particular, two-time multi-particle correlation and response functions for a one-dimensional 
diffusion-limited annihilation process with appropriate rates follow immediately from our results. Us- 
ing the exact mapping between diffusion-annihilation and diffusion-coagulation processes |23|. two- 
time solutions also follow for the latter processes. We are currently exploring this link, with the aim 
of e.g. obtaining exact scaling formulas for two-time correlation and response functions in the one- 
spin facilitated Fredrickson-Andersen model [24j at low temperatures. These should help to clarify 
the meaning and extent of the apparent trivial equilibrium fluctuation-dissipation behaviour found 
in numerical simulations |25L 
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APPENDIX A: PROJECTION OF THE GREEN'S FUNCTION 

Here we derive the identity 

E oB(t) 6 « = £ (-ir^ri^c^Ggs^w (ad 

that is used in Sec. Ill Bl for the inductive proof of (|32|l . The functions G^ k \t) and H n (t) are given 
in 130JI and IjClJI . respectively, and i <E N(k), J £ N(k — 2) are ordered indices. The matrix b^ k \ 
accounting for links between levels k and k — 2 in the hierarchy, is discussed in the text below (|1 T|) 
and may be written as 

k-i 

Now we focus on (| A 1 j) : for any fixed j' the summation variable j on the l.h.s. of (|A1|) must assume 
the particular values j = j_ £ in order to hit a nonzero matrix element of b^, where 

l<V<k-l: j n4 =(j' 1 ,...f n _ 1 ,£,£+l,j' 11 ,...j' k _ 2 ) 
V = k: j kjN = (l,j[,...j' k _ 2 ,N) 

The constraint j„ e S iV(fc) restricts the range of the parameters rj and I. If, for instance, 2 < r\ < k — 2 
the conditions j'_ x < £ and £ + 1 < j' must be satisfied. So only if j'_ 1 < f — 3 is there enough 
"space" to insert the indices {1,1+ 1) and retain strict ordering. Otherwise the corresponding vector 
j„t must be omitted. Nevertheless we may write 

.. j'l-i fc-2 i^-i JV-l 

i E 0)© = E G £U) + E E G S_w + G S fc , N w- ( A4 ) 

'ieJV(fc) i=i v =2e=f 7i _ 1 e=j' k _ 2 

Note that in contrast to what we have just said we insert the indices {£, £+1) between jL_x and j' in 
(j A4|) even if jL_x — jL — 1. This is possible because the Green's function £?W from 1)30(1 is identically 
zero if j n j contains an index pair, as follows from its permutational antisymmetry discussed after 
<|30|) . We may use permutational antisymmetry of G^ k ' once more to replace in (|A4Jl all j e with 
1 < r] < k by ji i since the corresponding permutations involve an even number of transpositions. 
The 77-dependencies in (| A4|) via Green's functions then drop out and we can combine the sums 

N-l N 

E Gg(*)6?J< = 7 E d (*) +7Gi3, w (*) =7E G S M w- ( A5 ) 

jeN(k) e=i e=i 

The second equality in (|A5|) follows when using permutational antisymmetry of and iV-(anti)- 
periodicity (|B9|) of the functions /„ (x) contained in G^"^ . The reasoning applied is in full analogy to 
the discussion of $( fe J in Sec. Ill Al Before proceeding, we discuss briefly why 1(A5J) is correct also for 
j'l = 1 or j'k-i = One of the sum over j 1 1 and j k _i 1 in l|A4(l is then empty and evaluates to zero, 
corresponding to the fact that j 1( G N(k) or j^-ie S AT(fc) cannot be achieved for any £ In l|A5jl. 
on the other hand, the corresponding terms drop out because j 1£ contains index pairs. 

The evaluation of IIABjl requires us to keep track of the first two components of j xt , viz. {£,£ + 1), 
in the Green's function (|30(1 . To do so we focus on the subset of permutations S^. u (k) C S(k) 
that map these first two components onto components /1, v 1 respectively. Any 7r € S^^k) may be 
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factorized into n = 7r Mjt/ off where 7f performs an arbitrary permutation on the components 3, . . . k and 
■Kfj, tV then puts the components 1, 2 into positions /x, z/, respectively. We have (—1)" = (— l) 7I " M ' y (— l) 7r , 
where (— 1) 7I >>" = (— \y-v- 1 for v > [i and (— l) 7r "'" = — (— for v < /i follows by counting the 
number of transpositions in ~K^v Using this factorisation for 7r G S fJi ^(k) and the fact Uu^w S^. v {k) — 
S(k) allows us to rewrite the Green's function (|30|l in the form 

i</i<w<fc 

- £ (-lr-^e-^/^-^TtJ^-J-iC^G^j^W. (A6) 

l<!/<fi<fc 

Relabelling n <-> in the second line of (|A6|) and substitution into (|A5|) then gives 

N 

j£N(k) l<n<u<k 1=1 

(A7) 

Our derivation of the identity IjAljl is completed by using the convolution property l|B8(l of the 
functions I n (x) to evaluate the £-sum 

N 

ie- 2t J2[li ll -i^-i-i-Ii.-iIi ll -l-i] (lt)= 1C - 2t [^- v -i-^- v+ i] (2<yt) (A8) 
and identifying i|A8() as the i-derivative of H^-i (2t) given by IjCljl. 



APPENDIX B: PROPERTIES OF 7„(x) 

The functions J n (a;) appear in our analysis of the Glauber-Ising model when we solve (|27[1 for $ 
and invert the Fourier transforms. We define 

T »to = ^ E z mq+xcosq = ^ E cos ( n «) e * COS9 < ( B1 ) 

where we have to set Q = Q G on even levels k of the hierarchy and Q — Q on odd ones. In the 
thermodynamic limit N — > oo the sets Q e ,Qo given in Q251 I26|) become dense on [0, 27r]. Hence for 
any fixed n£Z and iCKwe have 

lim J n (a!) = f p- cosMe 1 ™ 9 - I„(ar), (B2) 

that is the I n (x) on all levels k become modified Bessel functions I ra (x) as N — ► oo. For a comprehen- 
sive discussion of the I„ (x) see [2(j . Many properties of the I„ (x) in fact have an analogue for finite 
N. The following identities hold for Q = Q e as well as Q — Q : 

V0<n<N: I n (0) = 5 nfi (B3) 
VieK VneZ: J_„(x) = J„(x) (B4) 
Vx e R+ V - [f J < n < LtJ : ^0*0 > (B5) 

Viet VneZ: ^/„(x) = 1 [J n -i(x) + 4+^)] (B6) 
dx 2 

m+Af 



VxeR Vm,neZ Vq e Q o/o : ^ e lfc « J^°„ (x) = c - m «+* cos 9 (B7) 

771 + TV 

Vx,yeM Vm,neZ: £ /^(ar) I e J° k {y) = F n '°{x + y) (B8) 



k— 771+1 
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In IB 71 IB8JI the e/o symbols indicate that the identities hold on even as well as odd levels but not 
between even and odd ones. The functions I n (x) associated with even and odd levels differ in the 
following properties (here m,jieZ and x € R): 

Odd Levels (Q = Q a ) Even Levels (Q = Q c ) 

I n+N (x) = I n (x) I n+N (x) = -I n (x) (B9) 

m+N m+N 

^2 h(x)=e x l J *(a)|<e a (BIO) 

k—m-\-l k— m+1 

The identities fEU EJ EH follow immediately from the definition of I n (x). JEzHEEHEH, on the 

other hand, are a consequence of the convolution property and 7V-(anti)periodicity of the discrete 
Fourier transforms (|23l I24|) . We have not managed to prove the seemingly trivial inequality l|B5|) in 
a rigorous way. Numerical evaluations of (|B 1|) . however, leave no doubt that (|B5|) should be true in 
general. Assuming that l|B5(l is true, (jBlOfl can be verified as follows: on odd levels, where I n (x) is TV- 
periodic, (|B5|) implies that |/ ra (^)| = In{x) for all n S Z and hence l)B10jl follows from (|B7J| by setting 
q = 0; this is allowed since € Q . On even levels, the iV-antiperodicity of the I n (x) is cancelled by 
the modulus in (|B10|) . Hence we may shift the summation range to the region — L-^p'-J < n < LtJ 
where I n (x) is positive and drop the modulus. Then, by substituting (IB 111 and exchanging sums, the 
bound [{BlOjl can be verified [l^. To illustrate the actual shape of the I n {x) we show some plots in 
Fig.H 




-20 20 40 60 80 100 120 -20 20 40 60 80 100 120 



n n 

FIG. 4: Plots of the scaled functions e~ x I n (x) for a finite system of size iV = 100 (dots) and the N — > 00 
limit e~ x I„(x) (lines). The (left) right figure shows the iV-(anti)periodic functions associated with (even) odd 
levels, respectively. The curves correspond to x = 10 (narrow peaks) and x = 100 (broader peaks). 



APPENDIX C: PROPERTIES OF H n (t) AND RELATED FUNCTIONS 



The functions H n (t) originate from the projection of the Green's function (|AL|> . which is an essential 
ingredient for our solution (|40J) . Appendix lAl motivates the definition 



H n (t) = lj^ dr e- T [J»_i(7t) - Wtt)] 



(CI) 



When substituting i|Bl|l for I n (x) the integral in (|C1|) can be solved and we obtain the alternative 
representation 



H n {t) = i ^sin(nq) - 



7 sin g 
— 7 cos q 



1 - e 



-t (1— 7 cos q) 



(C2) 



In l|C2(l one has to substitute Q = Q e / as appropriate. Like the I n (x), the H n (x) thus have slightly 
different properties at even and odd levels of the hierarchy. 
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Two obvious properties of H n (t) are H- n {t) = —H n (t) and H n (0) — 0. In the thermodynamic limit 
N — > oo uniform convergence of I n (x) on any finite interval x £ [0, a] and for fixed n G Z allows us 
to replace / n (x) in (|C1|I by its limit l n {x). In (|C2|I . on the other hand, the limit N — > oo amounts to 
replacing by / dg. In analogy with Appendix 151 we use the notation H„(t) for the AT — > oo 

limit of H n (t). Further properties of the functions H n (t) may be obtained directly from 1|C1(1 using the 
features of I n (x) given in Appendix iBl It is interesting to note that, for 7 = 1, H„(2t) is essentially 
the probability that two independent, continuous time random walkers on a one-dimensional lattice, 
which are a distance n > apart at time t = 0, meet before time t 15]. 

In the context of Sec. lVUVll the following sum, which we use to extend the definition of H n to two 
time arguments, is of relevance 

Hfyh,t 2 ) =E e " tl ^n(7*i)<° n (*2). (C3) 

n 

Here e/o emphasises that both I and H on the r.h.s. of i|C3(l must be associated either with even or 
odd levels. From l|B8|) and IjClfl we immediately obtain 

tl+*2 

H n {tx,t 2 ) = \ J dre- r [^-i(7r)-/„+i(7T)] (C4) 
ti 

and thence, or from (|B7(I and l|C2|l . 

HJt u t 2 ) = 4 Vsin(ng) 7Smg e -«i (i-7cos g ) N _ e -t a (i- 7 cc«,)l _ (C5) 
AT 1 — 7 cos q L J 

According to e.g. 1C1IC4P we have the obvious link H n (t) = H n (0, t) between functions with one and 
two time arguments. But it is quite remarkable that we may conversely write 

H n (ti,t 2 ) = H n (h + t 2 ) - H n (h), (C6) 

when considering the definition l|C3J) . The expressions l|C4IC5|l are particularly useful in the thermo- 
dynamic limit N — > oo where (IC3(I becomes an infinite summation. In l|C4IC5J) . however, this limit 
is straightforward. As for (|C1I(J2|) we then use the symbol H„(£i,t 2 )- 

In Sec. IIVI we obtained another function H^(t) as given by (|52() that is closely related to H n (t). 
The physical relevance of 7J" (*) i s discussed in the text below (|52|l . In mathematical terms we may 
simply consider H'^(t) as a generalisation of H n (t) since we have H%(t) = H n (t) for 7; = 0. For 
7i 7^ 0, a representation of H'^(t) in the form IjClf) is not possible. Note that while H n (t) is defined 
on even as well as odd levels of the hierarchy, H^(t) is always associated with even levels, i.e. the 
sum in (|52|l runs over Q e . In analogy with (|C3|I we now extend the definition of H"(t) to two time 
arguments. The following sums are relevant for Appendix iDl 

HU(h,h) = Y, e ~ txI i-n{-rh)H'j_ n {t 2 ) (C7) 

H"^(ti,t 2 ) = ^sgn(j-n)e-* 1 7?_ n (7*i)^"-„(*2) (C8) 

Tl 

where 1 < i, j < N and the sum runs over n = 1, . . . N for a finite system. The superscripts e and o 
of I n (x) indicate that these functions must be associated with even and odd levels of the hierarchy, 
respectively. The sum <|C7t may be evaluated easily using l")2t and (|B7|I 



N 

_7i_sinq_ e _ t2 (i_ 7COS(Z ) 7 sing 
1 — 7i cos q 1 — 7 cos q 



1-e 



-t 2 (l— jcosq) 



(C9) 



In analogy to H we obtain from |K2"|l. fC^I and Jdjl the link = H'^O.t) and the non-trivial 

identity 

K{t u t 2 ) = K(h+t 2 ) - H n (h) (C10) 
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for the sum (fTTfj) . For 7; = the equations (fl?jl . (fC9|) and l(CTo| reduce to $£2&, J^U and 

(|C6ll . respectively, with i/" replaced by H. Consistent with our notation elsewhere we will denote 
the N -> 00 limit of ff" by H". 

In contrast to 1C7|I , the sum (|C8|) cannot be evaluated explicitly. But it may be rewritten in terms 
of the integral representation 

#;'(ti, <2) = sgn(n) fc) +7 / 1 dre- T /°( 7 T)ij;'(i 1 - r,t 2 ). (Cll) 

Jo 

Note that, in contrast to H'^, the function is even in n. As for (|CJT|) . (|C9|1 we will drop the 
double-primes in (|C8|) . IjClip for 7^ = 0, and denote the A — > 00 limit of H"(t 1 ,t 2 ) by H"(ti,t 2 ). D ue 
to uniform convergence of the integrand, this limit is obtained by replacing I,H with I,H in llClljl . 

Equivalence of (|C8Jl and (|C11|I may be verified as follows: substitute the expression for H"(t) 
into l|C8|) and exchange the sums over n and q. Now, the sum over n is of the form 

N 

fi,j(x, q) = s S n 0' ~ n ) T i-n( x ) sin [q (j ~ n)] , (C12) 

n=l 

where q is the summation variable in H" and thus q £ Q c - Differentiating (|C12|I w.r.t. x using (|B6|) . 
shifting the summation variable n — ► n ± 1 as appropriate and rewriting the result in terms of sums 
over the range n = 1, . . . A yields 

d 1 N 

0^ fiA*> l) = 2 £ s § n ^ ~ n + !) ^-n(^) sin[q (j - n + 1)] (C13) 
+ 5 s S n C? " iV ) ^-JV-iW sinfe (j - iV)] - i sgn(j) J?_ 1 (a:) sm(qj) 



1 * 

2 X! s s n ti - n - 1) ^-nW sin [<7 C? - « • - !)] 



n=l 

+ \ sgn(j - 1) J?(x) sm[q (j - 1)} - \ sgn(j - A - 1) /^(x) sin[g ( j - A - 1)]. 

Next focus on the second line in 1C13|) . Due to A-periodicity (|B9(I of the functions I n (x) on odd 
levels we have I°_ N _ 1 (x) — I°_ 1 (x). Also, since q £ Q c as given by J2EJ), the identity sm[g (j — A)] = 
— sin(gj) holds. Recalling, finally, that 1 < j < A is assumed and noting that sgn(j — A) = — sgn(j) 
for 1 < j < N thus shows that the two terms in the second line of (|C13(1 cancel each other for 
1 < j < N. When j = N both terms are zero because sin(0) = sin(Aq) = 0. It follows similarly that 
the fourth line in (|C13jl drops out. Expressing sgn(j — nil) via sgn(j — n) by adding compensating 
terms for j = n, n ± 1 and using the trigonometric identity sin(a + b) + sin(a — b) = 2 sin(a) cos(&) 
allows us to combine the sums in (|C13|I in the form l|C12|) . Consequently the fi,j(x,q) satisfy the 
inhomogeneous differential equations 

d 

fa fi,j( x > q) ~ cos 0?) fi,i( x > q) = fi-jW sin(g). (C14) 

From l|B3|) and IJC12I) we also have fi,j(0, q) = sgn(j — i) sin[q (j — i)]. Thus we obtain fij(x, q) by 
integrating i|C14(l for this initial condition, i.e. 

f i , j (x,q)=sgn(j-i) S in[q(j-i)}e xcos « + smq f Ax' e^') COB « J? , ,(x'). (C15) 



Using (|C15|I in (|C8(I immediately produces (|C11(I if one bears in mind (|C9(I . 

We conclude this section by showing that the general expressions l|C9(l , IjClljl simplify quite signif- 
icantly if one considers a quench from the random initial configuration corresponding to 7] = 00 to 
T = 0, i.e. 7; = and 7 = 1. Let us first focus on H^(t±,t2) as given by i|C9(l . which for 71 = we 
write as H n (t\,t2)- Since also 7=1, we may use the trigonometric identity 



sm.T , , . sim 



\(n + 1) x] = sin(?ix) h cos(nx) -I- cos[(n + 1) xl (C16) 

1 — COS X 1 — COS X 
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in (|C5(I . This yields the zero temperature recursion formula 

H n+1 (h, t 2 ) = H n (t 1 ,t 2 ) + e"' 1 [/„ + J n+1 ] (ti) - e -( tl+ *^ [/„ + 7„+i] (t x + t 2 ) (C17) 

for H n (ti,t%). Starting from Ho(ti,t 2 ) = 0, any H n (ti,t 2 ) can thus be expressed purely in terms of 
functions I n {x). By setting t% = in (|C17|) and using l|B3|) we also obtain a recursion for H n (t) = 
H n (Q, t), that is 

H 1 (t) = l-e- t [I +h](t) and tf n+1 (i) = H n {t) - e~* % + I n +i] (t), (C18) 

where 1 < n + 1 < N. Outside that range we have H- n (t) = —H n (t) and iV-(anti)periodicity in n 
according to e.g. (|C1|) . Equations (|C17IC18|) apply in the thermodynamic limit, too, if we replace H 
by H and I n {x) by modified Bessel functions \ n {x). 

Now we consider Hn(ti,t 2 ) at 7; — 0, 7 = 1 and, for the sake of simplicity, also N — + 00. As 
explained above, we drop the double-primes in (jCJl 1|) for 7; = and replace H,I by H,I in the 
thermodynamic limit. This gives 

iln(h,t 2 ) = sgn(ra)H n (ti,t 2 )+ f dr e~ T I„(r) H x (ti - r, t 2 ). (C19) 

Jo 

Expressing Hi(ti — t, t%) via (|C17|I allows us to rewrite l|C19|) as 

H„(ti,t 2 ) = sgn(n)H„(t 1 ,i 2 )+e^ 1 /„(t 1 )-e-(* 1+ *^/„(t 1 +t 2 ) 

+ e -< tl+t »> / drUtx + ta - t) [Iq + Ii] (r), (C20) 
Jo 

where the functions f n (t) are given by the single-sided convolution integral 

/„(*)=/ drl„(r) po + Ii](t-r). (C21) 
Jo 

By Laplace transforming l|C21() it is easy to verify that fo(t) — e*— Io(t) and f n +i(t) — f n (t) — 2l n+ i(t). 
This recursion formula produces a very similar expression for e~ t f n (t) as l|C18|) for H n (t) which in 
fact amounts to the link 

e-* l n (t) + e"< f n (t) = 5 n , + sgn(n) H n (t). (C22) 
Substituting l|C22l) into (|C20f) and using l|C6|l reduces the expression <jC19f) for H„(ti,t 2 ) to 

e- (l I„(^) + R n {ti,t 2 ) = c -(* 1+ * 2 ) |l n (i! + fa) +^ dTln(ti + * 2 - r) [I + Ii] (r)| . (C23) 



APPENDIX D: THE SUMS £ AND T 



Here we focus on the sums £, T given in l|61|). (|62|) . These appear in the expressions for the 
homogeneous solutions of the two-time correlation and response functions discussed in Sec. IVII 
Our intention in this section is to reduce the sums in 1|) . I|62|l over ranges scaling with TV to sums 
over a finite number of terms, which therefore remain manageable even for N — > 00. 

Let us first consider the one-dimensional sum £ given in (|61|l . where we use the abbreviation 
p = dim(_7'). For a finite system we have j € N(p) and the summation range is m = 1, 2, . . . N. If p is 
even, then the product over the sgn's in I jfiljl is +1 for m < j\ and m > j p . Therefore the expression 

dim(j) 

Sj,m = 1 ]J S S n (^ _ m ) ( D1 ) 
A=l 

can be non-zero only for ji < m < j p . In fact we have Sjj^ — 1 for v = 1, . . .p and Sj_ rn = 2 for 
J2r/-i < m < j 2v with 1/ = 1, . . .p/2. Otherwise 5j, m is zero. In terms of (|D1(I we may rewrite Ijfiljl 

as 

£ U (*- = E I 1 - e_A ^-™(7A) Hl_ m (2t'). (D2) 

m 
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Now recall that C^ k ' l \t,t') and R^ k ' L \t,t') are non-vanishing only if k + I is even. Hence if p = I (or 
p = I ± 2 for responses where from (|68[) the argument j in (|D2fl might be a vector j u ' s ) is even, then 
so must k be. Consequently the functions I n (x) in (|D2|> are associated with even levels. This allows 
us to rewrite l|D2(l 

(t, = (A, 20 - £ 5 jtm c- A / iE _ m ( 7 A) fl£_ m (2t') (D3) 

rn 

using (|C7|) . The remaining summation in (|D3|) only produces non-zero terms in the N- independent 
range ji < m < j p . Consequently the thermodynamic limit N — > oo of l|D3(l is straightforward. 

Following a similar strategy we now rewrite the sum £ for p = dim(_7') odd. In this case the product 
over the sgn's in gives +1 for m < ji but —1 for m > j p . Hence we introduce 

dim(j) 

<rj,m = s § n ( a -m)- Yl sga(jx - m), (D4) 

X=l 

where a is some reference site that should be chosen in the range ji < a < j p . Then (|D4|) again only 
produces non-zero values for ji < m < j p . In terms of (|D4J) and with the choice a = j„ the sum (|fj 1|) 
becomes 

C(*'*') = Eh n ^- m )-< 

m 

Applying the same reasoning as below (|D2J| we conclude that in (ID5|I the functions I n (x) must be 
associated with odd levels. Hence we utilise (|C8ll to rewrite (1D5|) as 

Sl dv (t, H) = Hl_ lE (A, 2f) - £ 4 m e- A 7„( 7 A) H'l_ m {2t'). (D6) 

m 

Because the summation range in (|D6|l that produces non-zero terms is again A-independent and 
bounded by j\ < rn < j p , the thermodynamic limit of l|D6(l is obvious. 

For the two-dimensional sum T given in (1621) with even p = dimj we may proceed in full analogy 
with £. Expressing the products over the sgn's in (|62|l by (|D1(I yields 

(*> = E t 1 - - + *i.*J e ~ 2A 7 ^-™(7A) /i,-«(7A) H';_ m {2t'). (D7) 

The unrestricted one and two-dimensional sums in (|D7|I may be evaluated using (|B8|I and l|C7|l . We 
are then left with bounded summation ranges ji < m,n < j p 

- 5 ^ C ~ A [^-m(7A) <_ m (A, 2i') - / is _ m ( 7 A) <_ ro (A, 2t')] 

+ E *i>™ e_2A 4— (TA) W7A) <'_ m (2f ). (D8) 

Rewriting the two-dimensional sum J 7 for p = dim(ji) odd is somewhat problematic. If we express 
the products over the sgn's in (|62|l using (|D4J| it is not possible to choose reference points a s.t. we 
arrive at sums of the form (|C8|) . It turns out that the most compact - yet still bulky - result follows 
when choosing 

T L ,i s (*> f ') = E [ s § n (^ ~ m) sgn(i s - n) - sgn(i e - rn) a? n - sgn(i 5 - n) erV m + crV m a? n 

x e- 2A / le _ m ( 7 A)4_„( 7 A)C m (20. (D9) 



e- A / 4e „ m ( 7 A)^_ m (2t'). (D5) 
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By generalising l|C12j) integral representations for the sums in 1D9(I may be derived. One should find 

?UW) = E4m4% e_2A ^- m (7A)^-„(7AX_ ro (2*0 

m.n 



^4 >ro e- A Ji e _ m ( 7 A) / dre- T J ( 7 r) [<_ m _x (A-r,2t') 
J° 

E 4%. C ~ A 4-™(7A) / dre- / (7t) [K-m-i - K-m+i] (A - r, 2f ) 
_ Jo 



7 2 f A - r 

— dne ri / (7n) / dr 2 e T2 / (7T"2) 
4 Jo Jo 

x [-^-i.- 2 + 2^_ Je - K-U+2] (2A - n - r 2) 2f ). 



(D10) 



As the reference points in (|D9|I are « £ and is the bounds on the summation ranges in l|PT0)l are of the 
form min(i e ,j'x) < m < max(i e , ji) etc. For 7; = and 7 = 1 the expression IjDlOfl may be simplified 
using (JOnjl. 



APPENDIX E: THE NOTATION OF EQUATION fTTb 

In this Appendix we illustrate by means of two simple examples how arbitrary homogeneous solu- 
tions for two-time multispin response functions may be obtained from (|71(l . 

Let us first consider the simplest possible case of the spin response function to a local magnetic 
field. Here we have k = I = 1 and i = (i), j = (j). According to the formalism presented in Sec. IVII 
the first step in the evaluation of l|71|) is to determine the squashed index vectors Ijfi8(l ■ In the case at 
hand there is only one such vector since I — 1, viz. j ' s = j 1 = (j — + 1) with Z 1:S = 3. Then, 
by setting k = I = 1 in l|71|) and substituting I 1 ' 3 , we obtain 



-) E (-irn v ' 

' tt67"(1) A=l 



(«Uj')„(2A_l),(iUj)„(2X) 



(t,0 



Y E ( 1 )' r II 1; 7<Uj 1 ' s )„(2>_i),(*Ui I '") w(3J . ) (*' t ') 



7rSP'(2) 



The summations over pairings in this equation are restricted to those, V, that contain the pair 
where we set V[ 7 - = 1. In the first sum we have 'P(l) = {Id} and pairs are drawn from iU j = 
Hence T"(l) = {Id} and the only term in this sum is V,' ■ = 1. In the second sum, however, we have 
to consider the 3 pairings (1,2,3,4), (1,3,2,4) and (1,4,2,3) contained in V{2). As pairs are drawn 
from i U ji 1,s = — + 1) the only pairing contained in V'(2) is (1,3,2,4); the sign of the 

corresponding permutation w is (— l) ff = — 1. So there is again only one term in this sum, namely 



"Vi- 



1,3+1) 



and therefore 



(1,1) 



( 7 A)[(l-^) 



¥)-^v/:; j+1 (m') 



Since both indices of V are components of j ' s we have to use the (a, b) = {j^jv) case of H6()|) to 
express this function. Making the corresponding substitution in the above expression produces the 
result G3 given in Sec. IVII 

The next higher nontrivial response functions are those with k + 1 = 4. As a second example we 
consider the case k = I = 2 where i — {11,12) and j — We restrict ourselves further to 

vectors j that satisfy ji = ji + lj corresponding to a perturbation that couples to adjacent sites of 
the lattice. In this situation we have two index vectors j v given by j = (ji — l,ji,ji + 1, J2) and 
j 2 = {ji,]2 — 1, J2,j2 + 1)- From (|S5|) we obtain the corresponding squashed versions j ,s = (ji — l,ii) 
and j ,s = (j2,j2 + 1) using j'2 = ji + 1. Their dimensions are obviously i 1,s = l 2,s = 2. Now, by 
substituting k — I = 2 in l|71|l. explicitly writing out the summations over \i and v and plugging in 



34 



l v,s we obtain 

$ (2,2) 



( 7 A) 



tt£P'(2) 



A=l 



y e (~ i ) 7r n v 'oiuj i - 3 )„ (2A _ 1) ,(iuj i . B ), (2A) 

7r£7"(2) 



(t,0 



A=l 



( 7 A) 



1_ t) 51 ^ _1 ) ,r II V '( : iUi), (2X _ 1) ,(iUi) T(2A) 
7r£7"(2) 



(t,0 



A=l 



"T E ( 1 )"II l; '{iUj 1 ' s ),( M - 1) ,(iU; 1 ') I ( 2 x)(*' t ') 



ttGP'(2) 



e- A / n _ 32 ( 7 A) 



A=l 



1-t) E (-^Ilnv 

7rGT"(2) 



(2X-1) 



>(*ui), 



(2A) 



(t,0 



A=l 



ttGP'(2) 



A=l 



e- A / i2 _,- 2 ( 7 A) 



1_ t) X! (~ 1 ) 7T il V '( J iUj), (2A _ 1) ,(iUj)x(2A) 

7r£7"(2) 



(t,0 



2 E nrn^^ 

7rG7"(2) 



A=l 



)„(2A-l),(iUj 2 ' H ), 



,(M') 



A=l 



Within the 4 square brackets the pairings are restricted to those, V, that contain the pairs ji), 
(*2, ji), (iiiia) an d {^2,32) from top to bottom, respectively. The functions V' carrying these indices 
are to be replaced by 1. In all pairing sums the sets V'(2) are of course subsets of V{2)\ the latter 
contains the pairings (1,2,3,4), (1,3,2,4) and (1,4,2,3). We now focus on the first square bracket; 
the remaining ones follow similarly. In the first sum, drawing pairs from i\Jj = (ii, 12, ji, J2), the only 
pairing containing is (1,3,2,4) and has the sign —1. So this sum reduces to — H*i ji^i 2 ,i 2 = 



—Vi 2j j 2 . The second sum, on the other hand, draws pairs from iUj :S = (ii,i2,ji — 1, Ji) ■ Hence the 
only pairing that contains is (1,4,2,3) having the sign +1. This sum therefore collapses to 

+1^ 31 K' 2 = +Vi 2 .j 1 _i. Applying this reasoning to all summations yields 



*( 2 : 2 > = e - A I„ 



+ e- A L, 



-A 



1 - 



v', 2 (M0- VCU(M') 



i-~4 v. 



1, 



12- J2 \ 



(7 A) -1-L W,(M') + 1V ! 



2 K ii,j 2 + 1 



(f,f) 



Note that the indices of all functions V comprise components of i and j or i and j v,s . So all of them 
must be expressed via the (a, b) — {i e ,jv) case of (IrJUjl . We remind the reader that for in (— l)"^ 1 
in H6()(l one has to substitute the actual component number of the index j v — (j)„ in j. This means, 
for instance, ji — 1 = (j' 1,s )i is component 1^ = 1 in j l s while j\ = [j ' s )% is component v — 2 in j ,s . 
Accurately applying (|60|l to the above expression yields the result l|74|) given in Sec. I VII 
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